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Phase transition from hadronic matter to quark-gluon matter is discussed for various regimes of 
temperature and baryon number density. For small and medium densities, the phase transition is 
accurately described in the framework of the Field Correlation Method, whereas at high density 
predictions are less certain and leave room for the phenomenological models. We study formation 
of multiquark states (MQS) at zero temperature and high density. Relevant MQS components of 
the nuclear matter can be described using a previously developed formalism of the quark compound 
bags (QCB). 

Partial-wave analysis of nucleon-nucleon scattering indicates the existence of 6QS which mani- 
fest themselves as poles of P-matrix. In the framework of the QCB model, we formulate a self- 
consistent system of coupled equations for the nucleon and 6QS propagators in nuclear matter 
and the G-matrix. The approach provides a link between high-density nuclear matter with the 
MQS components and the cumulative effect observed in reactions on the nuclei, which requires the 
admixture of MQS in the wave functions of nuclei kinematically. 

6QS determine the natural scale of the density for a possible phase transition into the MQS 
phase of nuclear matter. Such a phase transition can lead to dynamic instability of newly born 
protoneutron stars and dramatically affect the dynamics of supernovae. Numerical simulations 
show that the phase transition may be a good remedy for the triggering supernova explosions in the 
spherically symmetric supernova models. A specific signature of the phase transition is an additional 
neutrino peak in the neutrino light curve. For a Galactic core-collapse supernova, such a peak could 
be resolved by the present neutrino detectors. The possibility of extracting the parameters of the 
phase of transition from observation of the neutrino signal is discussed also. 
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I. INTRODUCTION 

The problem of transition from the confined phase of hadronic matter to the deconfined quark-gluon phase is 
now widely discussed both on theoretical [l| and experimental levels [H and also using lattice simulations Q. The 
transition process is nonperturbative in nature and requires for a realistic description the use of nonperturbative 
models. Nambu-Jona-Lasinio (NJL) model and its generalizations describe successfully the spontaneous breaking of 
chiral symmetry but, unfortunately, do not lead to confinement and can not tell much about the deconfinement phase 
transition. 

Recently the Field Correlator Method (FCM) Q , successfully applied earlier to the description of confinement and 
hadrons in the vacuum, has been extended for the description of hadron matter at nonzero temperature and density 
including phase transition into the deconfined phase. In this method, the quark-hadron phase transition occurs due to 
in-medium modifications of the QCD vacuum. The energy density of the QCD vacuum is determined by the gluonic 
condensate, and the colorelectric component of it is responsible for the confinement. In course of the quark-hadron 
phase transition, the colorelectric component evaporates, whereas the colormagnetic component stays intact (for a 
review see [5|). 

In this way one obtains the critical transition temperature as a function of the quark flavors in agrement with the 
lattice data Q. The phase diagram extrapolated to small temperatures points toward the critical quark chemical 
potentials of order /icr — 0.6 GeV. This value is large enough to allow the ordinary nuclei to exist as collection of 
nucleons rather than quarks. The low-temperature phase transition gives rise to strong nonperturbative attraction 
in colorless channels [3] that provides dominance of the qq correlations over the diguark qq ones, making thereby 
conjecture on color superconductivity not quite realistic, in line with earlier findings 8]. In two-color QCD, where qq 
and qq are equivalent, sup erconducting phase of quark matter with colorless quark Cooper pairs [9l4ll| is confirmed 
by lattice simulations 

The low-temperature region of the phase diagram is, however, less certain because no interaction in the hadron 
medium (nuclear matter) was taken into account yet. The standard models of nuclear matter (see e.g. jisj ) based on 
the realistic NN potentials predict energy per nucleon E/A depending rather sensitively on details of NN interaction. 

An additional uncertainty comes from the forces and nucleon-hyperon transitions, which make computation of 
equation of state (EoS) of the nuclear matter a rather complicated task, which calls for experimental verification. At 
this point an important piece of information can be obtained from astrophysics. 

During the last decade, great progress is made in observational astrophysics towards the study of properties of 
neutron stars (see e.g. and, for a review, EoS of nuclear matter determined by the density dependence 

of E/A is important for calculation of structure and cooling rates of neutron stars and possible appearance of quark 
matter in interiors of neutron stars as was suggested for the first time by Ivanenko and Kurdgelaidze [l7j . The current 
status of the quark stars hypothesis is discussed in Refs. [l6l.lT8|. 

The recent interest to the nuclear matter EoS is also connected to new constraints obtained from the collective flow 
data and subthreshold kaon production in heavy-ion collisions [l^, H^l . 

One of the main purposes of the present paper is to consider a generalized picture of TV TV interaction and of nuclear 
matter, which includes from the very beginning the quark degrees of freedom in the form of MQS. This implies a 
revision of the standard picture, based on the Yukawa mechanism of meson exchanges. 

In 1935 Yukawa [llj proposed a hypothesis that the interaction between the nucleons may be due to the exchange 
of finite-mass meson. The experimental searches have led to the discovery of pious, heavier mesons, and eventually 
to the emergence of One Boson Exchange (OBE) model of nucleon-nucleon interaction. In this model, pious are 
responsible for long-range part of the nucleon-nucleon potential, while heavy mesons are responsible for the interaction 
at intermediate and short distances ,2^] . 

OBE models provide an accurate quantitative description of nucleon-nucleon interaction. The physical meaning of 
the meson exchange at short distances is, however, not entirely clear because of the finite size of nucleons and mesons. 
The values of the proton and pion charge radii < >i/2= o.875 ± 0.007 fm and < >i/2= o.659 ± 0.025 fm ^ 
indicate that the mechanism of pion exchange is justified at distances r > 1 fm while at smaller distances there is 
overlap between hadrons, where the quark-gluon degrees freedom come into play and affect the dynamics. 

Restrictions of this kind should be taken into account when calculating EoS of nuclear matter. The maximum 
density of nuclear matter below which the OBE models can be applied is determined by comparing the proton charge 
radius with the average distance between the nucleon with its nearest neighbor. The nearest neighbor is localized 
between r and r -I- dr, while the sphere of radius r with the probing nucleon in the center is empty. The probability 
distribution of the nearest neighbor is the probability of not finding a nucleon inside the sphere multiplied by the 
probability of finding a nucleon in the volume element dV — Airr'^dr. The probability do not find a nucleon inside 
the sphere is given by the Poisson law Pq — exp(— pV^), where p is nuclear matter density. For the saturation density 
po = 0.16 fm~^, a simple calculation gives < r >= 1.02 ± 0.37 fm, where < r > is the mean distance, the second 
number is the standard error. When the density is close to saturation density, OBE mechanism, obviously, can not 
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be applied. 

The problem of nuclear matter EoS is studied in the realistic OBE models of nucleon-nucleon interaction since the 
mid-1970's. One can distinguish three main approaches: A popular class of models based on the mean-field approx- 
imation [23, HI]. The field-theoretic Dirac-Brueckner-Hartree-Fock method goes beyond mean-field approximation 
(see and references therein). Variational method is described in Ref. [26| . 

One should note that in all approaches of OBE-type additional fine tuning is required to obtain the realistic values 
of energies and densities, e.g., in the variational approach, 3N forces have to be added to the NN potentials. 

OBE models are in reasonable agreement with laboratory data [l^dSlj but predict a surprisingly low value of the 
maximum mass of neutron stars in the /3-equilibrium, if one includes hyperons [13-1111 ■ /3-equilibrium leads to the 
occurrence of the hyperons when density increases to about (2 ~ 3)pq and, consequently, leads to a softening of the 
nuclear matter EoS. 

It has earlier been noted [HI that observational data on the rotation speed of X-ray transient XTE J1739-285, 
which point to a soft EoS, and mass of the pulsar PSR J1748-2021B, which points to a very stiff EoS, are almost 
mutually exclusive. In a recent study [s^] data on compact sources 4U 1608-52, 1820-30 and 4U EXO 1745-248 were 
re-analyzed. The authors came to conclusion that only a very soft EoS is consistent with values of the mass and 
radius of the stars. A similar technique was used previously in the analysis of data from a compact source EXO 
0748-676. In conjunction with the analysis [lH|, one can assume the existence of two classes of compact stars, e.g., 
neutron stars metastable against conversion to exotic stars, such as quark stars, strange stars, or perhaps dibaryon 
stars. The results of Refs. [15:, JJSj still require confirmation. 

The difficulty in describing the massive neutron (hyperon) stars as well as the possibility of existence of various 
classes of compact objects stimulate the search for new concepts and models in which quark-gluon dynamics at short 
distances plays a more prominent role. One can expect that the microscopic models are more adequate at high density 
and can provide a physically satisfactory picture of the short-distance dynamics and high-density EoS. 

The standard realistic NN forces (and also 3A^ forces) used for nuclear matter calculations exploit only baryonic 
and mesonic degrees of freedom. However, multiquark states (MQS) can appear in nuclear matter also, yielding 
degrees of freedom of their own. 

This idea is discussed, e.g., in Refs. [33 - l36| . In most papers, it is associated with multiquark bags similar to 
the MIT bags [s^l, where only perturbative interquark forces act inside bags. This invokes immediately the idea of 
deconfined pieces of matter inside of hadrons and the quark-hadron phase transition caused by the increased number 
of multiquark bags and their final overlap. 

However, the MIT bag is only a crude model of hadrons. Lattice and other analytical models show that the 
correlation length in the QCD vacuum is very small ~ 0.1 fm. This fact is due to large mass of nonperturbative 
structures called gluelumps (for a review see [38|). The string (confinement) between quarks acts therefore already 
at distances ~ 0.1 fm or less and instead of multiquark bags one gets for MQS strongly bound systems with inner 
density growing with the number of quarks and a radius twice smaller than in the MIT bag model [39| . Moreover, 
a detailed study of the dependence of the confinement potential on the density [13] showed that the medium effects 
lead to an additional attraction, as a result of which MQS are becoming smaller and lighter. A growing admixture of 
MQS can produce another minimum in the curve of E/A vs density providing a phase transition into the heterophase 
nuclear matter with substantial MQS component or the quark matter. 

MQS can play an important role in the so-called cumulative and subthreshold processes i.e., reactions on nuclei that 
cannot proceed on single nucleons 41|. The corresponding cumulative number Ncum defines the minimal number of 
nucleons needed kinematically for the reaction. It is clear, that MQS are appropriate objects to provide the cumulative 
effect. The MQS density enters directly cross-sections of the cumulative and subthreshold reactions. There is a close 
connection between the MQS density in nuclei and neutron star interiors. 

Multiquark configurations that appear first when overlap of nucleons becomes significant can be considered as a kind 
of dibaryons. The experimental searches of dibaryons in the past did not give conclusive results. Recently, resonance 
behavior of the double pionic fusion reaction pn — )■ dir^TT^ measured at CELSIUS- WASA has been interpreted as 
evidence for a A A dibaryon [43|. Such dibaryon has, however, more features in common with deuteron rather than a 
compact 6QS. 

The possibility for occurrence of a Bose condensate of dibaryons in nuclear matter is discussed in Refs. [sj, IH, |43| - 
|48| . The ground state of the heterophase nucleon-dibaryon matter is stable for a wide range of parameters of the 
models discussed and nuclear matter densities. Using the mean-field approximation of the OBE models, constraints for 
the w- and cr-meson coupling constants with dibaryons were extracted from properties of nuclear matter at saturation, 
stability condition of the binary mixture of nucleons and dibaryons, and from the existence of massive neutron stars. 

To give a quantitative method for study of the nucleon and dibaryon fields, we are using the so-called Quark 
Compound Bag (QCB) model [401, appeared in the development of the P-matrix formalism of Jaffe and Low |5C| . 
Jaffe and Low proposed to identified MQS with the so-called "primitives" which appear as poles of P matrix rather 
than S matrix. The P-matrix formalism was used to yield an accurate description of the nucleon-nucleon systems in 
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Ref. [51| (see also (49|, l52|-l55|). 

The dynamical character of the QCB model allows to apply it for the description of nuclcon-nucleon interaction at 
finite density and temperature and for the study of nuclear matter EoS. 

It should be noted also that the QCB model is very economical since it uses only a few parameters, such as mass 
and radius of the MQS with a fixed orbital and total angular momentum of quarks. 

The average orbital momentum of a nucleon pair in the nuclear matter can be estimated as L < rpp — 1-4 ± 0.5. 
This estimate does not depend on the density. The s-channel dibaryon exchange as the mechanism of nucleon-nucleon 
interaction is therefore restricted in nucleon matter by a few lowest partial waves. An additional suppression of the 
high partial waves appears in the low-density limit because of the threshold behavior ~ a{pF)^ of the partial wave 
amplitudes, where a is the scattering length. The S'-wave interaction is therefore dominant in the low-density limit. 

The comparison of the observational data with the calculations of static properties of neutron stars and the dynamic 
simulations of core-collapse supernovae rise questions which have not found a definite answer yet. One such issue is 
the quantitative properties of the phase transition. The inclusion of MQS in the NN interaction dynamics, e.g., in 
the QCB model is an alternative approach to describe nuclear forces and nuclear matter EoS. In this approach, the 
modeling of heterophase states of high-density baryonic matter is possible. 

The physical conditions that appear at the final stages of the core collapse in interiors of newly born hot protoneutron 
stars are unique to constrain parameters of the quark-hadron phase transition, which cannot be constrained otherwise 
from laboratory data (see, e.g., (56j). Observation of neutrino bursts from supernovae, which may genetically be 
related to a phase transition, will play an important role in the identification of physical nature and properties of the 
exotic states of nuclear matter at high density. 

Massive stars end their life with gravitational collapse of iron central cores that according to astronomical obser- 
vations must result in the explosion of so-called core-collapse supernovae (SN). However, an extensive hydrodynamic 
modeling during already more than thirty years has demonstrated that in case of spherical symmetry it is very hard 
to simulate the explosion: stellar envelope would not separate and finally fell back on the collapsed core. Therefore, 
the SN theorists began to concentrate special attention on basically nonspherical effects such as rotation, magnetic 
fields, large scale convection, jet streams, etc (e. g., see [13, [13 and references therein). 

Here, we discuss the calculations of gravitational collapse that was undertaken to estimate the effect of the nuclear 
phase transition on the core collapse dynamics. Although the property of phase transition to destabilize the hydrostatic 
equilibrium of stars and giant planets is well known for a long time [59'-'63'| the hydrodynamic consequences of such 
a destabilization were yet rarely addressed in the investigation of SN mechanism. For a historical review and as the 
starting point for understanding the instability induced by phase transition, remarkable essays by Seidov 
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be recommended. It was understood that the phase transition onset in stellar center can generate a shock wave 
However the hydrodynamic effects on the SN mechanism expected from such a shock required further detailed study. 
Recently there appeared a detailed research of possible phase transition influence on the dynamics of the collapse 
and SN mechanism [68l [69| . The research is based on sophisticated hydrodynamic code including neutrino transport 
and quark-hadron phase transition simulated with the aid of the MIT bag model. It was shown that an additional 
neutrino peak appears in the neutrino light curve as a specific signature of the phase transition. For the core-collapse 
supernovae of our Galaxy, this peak could be resolved by the existing neutrino detectors. 

The QCD phase transitions may be a good remedy for triggering the SN explosion in spherically symmetric SN 
models, while the detected neutrinos can provide important information on physical nature and parameters of the 
QCD phase transition. 

This work reviews phenomenological and microscopic approaches to NN interactions and nuclear matter properties 
with inclusion of quark-gluon degrees of freedom, quark-hadron phase transitions and influence of phase transitions 
on physics of neutron stars, supernova outbursts and neutrino signals from supernovae. 

In Sect. II we discuss the — T QCD phase diagram, give estimates of critical temperature of the deconfinement 
phase transition at low baryon chemical potentials, and discuss the role of diquark correlations in cold quark matter. 
In Sect. Ill we discuss modification of the confinement force at nonzero chemical potential and estimate critical 
density above which formation of MQS becomes energetically preferable. In Sect. IV the QCB model is introduced 
in its original nonrclativistic form, while Sects. V-VII are devoted to the relativistic formulation of QCB and the 
accurate description of NN data in the QCB framework. In Sect. VIII the EoS of nuclear matter are given in the 
standard OBE form, in Sect. IX in a simplified nonrclativistic QCB form, while section X is devoted to the relativistic 
mean field treatment based on the QCB model of Sect. V. In Sect. X we provide equations to calculate the G-matrix 
and describe qualitative features of the phase transition to the new MQS phase. In Sect. XI astrophysical aspects of 
(proto) neutron star dynamics are given and the evolution of the supernova outbursts is presented based on detailed 
calculations. Section XII is devoted to conclusions and perspectives of the present approach. Appendices contain 
technical details necessary to derive some of the equations in the text. 
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FIG. 1: QCD phase diagram in the temperature-chemical potential plane. First-order phase transitions are shown by solid 
curves. Filled circles are second-order phase transitions. Crossover at /i « shows transition from hadronic phase to quark-gluon 
plasma. Crosses depict heavy-ion collision experiments. 



II. QUARK-HADRON PHASE TRANSITION IN QCD 



In Fig. [l]one can see the generic picture of the phase transition expected in QCD. Only the low-/i region is accessible 
to lattice QCD, while mostly NJL-type models are used to predict the phase curve and the superconducting phases. 
We shall show below that the phase curve can be obtained from the Field Correlator method (FCM) [1, 0| and indeed 
looks like in Fig. [H whereas the superconducting phases are unlikely. We shall also find the critical values Tc and 
for symmetric nuclear matter {p,p = fj,n = 3/ig) with switched off weak interactions in terms of the gluon condensate, 
following 0. 



A. Gibbs' criterion 



According to the Gibbs' criterion, two phases in thermodynamic equilibrium have the same temperatures, pres- 
sures, and balanced chemical potentials. Phase with highest pressure is thermodynamically preferred provided other 
conditions are equal. Chemical equilibrium leads to a set of relations for the chemical potentials of constituents. If 
weak interactions are switched off, chemical equilibrium for the substance with Uf flavors reduces to 

flA = C^fla, (II.l) 

where and are the chemical potentials in the hadron phase with A = p,n, A . . tt*- and Tr^-mesons, etc. and 
in the quark phase with a = u-, d-quarks for Uf = 2 and u-, d-, s-quarks for n/ = 3. The matrix determines the 
composition of particles in the hadron phase and, in particular, allows to express electric charge of particles in the 
hadron phase ca m terms of electric charge of the constituents Ca- ba = C\ea- If leptons appear in the substance 
then A and a run over the leptons too. For a representative set of hadrons and electrons and muons has the form 

a — )■ u d s e fi 

2 10 

1 2 

2 10 

-1 1 • 

1 

1 



A P 
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The critical temperature without the lepton component depends on the n/ chemical potentials of the quarks, that 
are independent parameters. 

If the relaxation time is longer than the typical time scale of the weak processes, the weak interactions come 
into play. This is the case of phase transitions in neutron stars whose lifetimes are comparable with the age of our 
Galaxy. Even in supernovae the free fall time is much longer than time needed to keep substance under the chemical 
equilibrium with respect to the weak interactions. 

The weak interactions modify the matching conditions (III. II) . The lepton component appears in the substance and 
the electric neutrality is to be imposed. The modified conditions look like (see e.g. [7C|) 

MA - C'X^la + caV, (II.2) 

where ca is electric charge of particle A, the indices A and a run over leptons also. The chemical equilibrium generates 
at the boundary of two phases a jump of the electrostatic potential V formed by a double layer of leptons to match 
chemical potentials under the condition of the bulk electroneutrality. The similar jump exists e.g. at boundaries of two 
different metals in contact. The chemical equilibrium with respect to the weak interactions and the electroneutrality 
condition leave each phase with one free parameter only that can be chosen to be, e.g., the neutron chemical potential 
in the hadron phase and the d-quark chemical potential in the quark phase. The matching conditions pi.2[) are 
reduced to one equation for and at the phase transition. The critical temperature is thereby a function of one 
parameter only. 

Below we start with the system of quarks and gluons without weak interactions and where interaction in the hadron 
phase is also switched off, as it happens in the large Nc limit. One can as a first step neglect also interaction between 
quarks and gluons as compared to their interaction with vacuum fields. This approximation can be called the Vacuum 
Dominance Model (VDM). 



B. QCD phase diagram in FCM 

The basic notion for the phase transition in QCD is the nonzero vacuum energy density Syac introduced in QCD 
sum rules j7i| . The conformal anomaly gives 

Svac - 1/4«MM - 15q; \y^t^i^> I = ^ '^2 

where [ll El 

G2 = (0.01 ± 0.002) GeV''. (II.4) 

The quark- hadron phase transition occurs due to reconstruction of the vacuum in the course of which the colorelectric 
part oisyac, ensuring confinement, is thermodynamically not advantageous and therefore vanishes in the temperature- 
affected QCD vacuum fz^. 

In the confining phase the pressure looks like 

Pi{T,fii) = e,ac + T\i{T,fii), (II.5) 

where denotes the set of the chemical potentials, the second term stands for pressure of the hadron gas. In the 
deconfined phase, the colorelectric confining correlator D^{x) — [73], as confirmed by lattice data [tJ, so that the 
vacuum energy density £,y^^ is decreased by about a factor of two, AG2 ~ within 10%. 
The pressure of the deconfined phase equals 

P2(T,fi2) = stl + T\2{T,fi2) (11.6) 

where 1^2 is the set of the quark chemical potentials, X2 (2^7/^2) = Pgi +Pq is a function weakly depending on T for low 
chemical potentials, pi — Pi/T* are reduced gluon and quark pressures. 
The Gibbs' criterion, 

Pl{T,,^Ji^) = P2{T,,^i2), (11.7) 

together with the matching conditions (jll.l[) allow to find Tc in terms of Ae„ac ^ AG'2. By neglecting the hadron 
pressure xi(r, /J,2), which is 0(10%) of Asuac for low chemical potentials, one gets 

Ae ^ 



T, = 7^""'- ^ . (II., 

\X2(Te,A^2) ' ^ 
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Xi(T, /Lt2) is, however, important for the order of quark-hadron phase transition. 

Next step is the evaluation of Pq and Pgi . It is shown in [ 76.] , that in VDM the main dynamical contribution to quark 
and gluon pressure can be written in the form of the Polyakov line (|IL9|) . creating effective selfenergies Vi{oo,T) for 
quarks and jVi{oo,T) for gluons, and otherwise quarks and gluons are free. The resulting expressions for Pgi,Pq are 

in (jll.lOl) . (jll.lip . Neglecting quark mass one obtains the form (jII.12[) . where 'I'^'L are given in (IIL13p . 

Equations (jll.lOp . (|II.12[) contain all necessary dynamical information for quark-hadron phase transition and dy- 
namics of quark-gluon plasma (QGP). Basic quantity is the potential Vi(oo,r), entering Polyakov line (jIL9|) . which 
is calculated on lattice and analytically see 



\L fundi = exp 



'2T 



(11.9) 



16 , , 16 
Pgi = —LadjiT) = — exp 



8T 



(11.10) 



7r2 ^ n4 

n—l 



fj,n 



(11.11) 



Pq{mq = 0) = 



12n/ ^ (-) 



n+l 



.i»cosh^ = 



7^2 „4 / y 

n=l 



(11.12) 



X dx 



(z) 



X dx 

oX + Z I 



(11.13) 



Now one can compute Tc(/i) for different Uf from Eq. (|II.8p . where ^AG2 
In this way one obtains a simple formula for the transition temperature 



(11-1 



32 



Go and 



Pq ~ Uf 



(|IlTT1) 



(11 



(11.14) 



where Pgiipq) are gluon (quark) pressure divided by in the deconfined phase, calculated in the same method (see 
for details). This simple formula yields for different number of flavors Uf correct values Tc{nf), in good agrement 
with lattice data f6|. 

In doing so one assumes that temperature and moderate density do not affect strongly the vacuum fields in G2 and 
for the latter one can use the same phenomenological values as for zero temperature, i.e. G2 = 0.01 GeV*. 

Proceeding in this way one obtains the phase diagram in the ^q — T plane @, see Fig. [31 where the critical value 
Mcr = Mg(^ = 0) — 0-6 GeV. This value is large enough and it will be shown below that vacuum structure is affected 
by the chemical potential /i^ and confinement is modified. As it is, one has for small T the first order density phase 
transition to the deconfined quark-gluon matter with strong np attraction in the white systems [3] preferring qq system 
over diquarks and thus making the proposal of quark superconductivity not very realistic, in line with earlier criticism 
of [8^. 

Then one obtains the numbers in Table HI where the asterisk corresponds roughly to 1/2 of standard condensate. 
The values Tc{nf = 0) = 0.27 GeV and Tc{nf = 2) = 0.19 GeV are in good agreement with numerous lattice data, 
while /ic w 0.6 GeV cannot be obtained on the lattice. The resulting curves Tc{ij,) are shown in Fig. [21 

The critical density of the phase transition into the strange quark matter (n f — 3) at zero temperature gets positive. 
This means that the strange quark matter is not absolutely stable. The condition of absolute stability of the strange 
quark matter restricts the critical temperature for vanishing chemical potentials by Tc{nf = 3) < 122 ± 7 MeV [Tq . 
According to the lattice data, this restriction is not satisfied with a large certainty. 
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TABLE I: Critical values of T and ^ for various n / and AG2 

★ 



AG2 

0.01 GcV"* 


0.191 


0.341 


0.57 


1 


rc( GeV) n/ = 


0.246 


0.273 


0.298 


0.328 


rc( GeV) n/ = 2 


0.168 


0.19 


0.21 


0.236 


T,{ GeV) n/ = 3 


0.154 


0.172 


0.191 


0.214 


Ale ( GeV) n/ = 2 


0.576 


0.626 


0.68 


0.742 


/ic( GeV) 71/ = 3 


0.539 


0.581 


0.629 


0.686 



In the derivation of the phase curve in Fig. [5] it was, assumed, that the pressure of the hadronic phase can be given 
with good accuracy by the vacuum fields only (the Vacuum Dominance Model (VDM)), while the contribution of the 
hadronic gas is neglected. This approximation is reasonable for /i = 0, where the effect of pionic gas is around 10%, 
but maybe crude for large /x, where hadronic matter is dense nuclear matter. Therefore the region of large /i and small 
T can be modified, when interaction in nuclear matter is properly taken into account. Since the main emphasis of 
the present paper is the treatment of the possible new phase of MQS at smaller /i, it will be reasonable to reconsider 
the hadron-quark phase transition with account of possible MQS contribution. 



0.2 

0.15 



0.1 



0.05 




0.1 0.2 0.3 0.4 0.5 0-6^ 



FIG. 2: Phase transition curves 
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C. Suppression of color superconductivity 

As was shown in [tJ [iH the diquark nonperturbative interaction is expressed via Vi (r, T) as follows (for a large 
size superconducting diquark) 

VQQ{r) = ^^1(00, T) + ^Vi{r,T) -> Vii^.T,) « 0.5GeV (11.15) 



For each quark it gives the factor 



^,m-exp(-^%^). (11.16) 



Hence each quark carries Boltzniann factor 



L, « exp ( ) . (11.17) 



Assuming that the dependence on /i is weak one expects that approximately the same factor holds for fi — > fXcrit- 
For the white QQ system this factor is missing [t^, hence one expects 

^.i».„p (-_)«! („,8) 

Therefore one expects QQ pairing, and not QQ pairing, in qgp. 



D. Critical baryon chemical potential in cold nuclear matter 

One can expand (|II.13|) for T 0, a — 00 



For the pressure one obtains 



and for the density 



+ + + a ^ (11.19) 

-V ''4 2 60 ' T ' ^ ' 



V 7r2 1 V 2 J V 2 
The critical /x at small T has an expansion 



^.(T ^ 0) = ^ + (48)V4rW 1 - ^- L-—+ . (11.22) 



2 (Mc - ^) 



For Vi{Tc) = 0.5 GeV and n/ ^ 2 one h as: ^lc{T = 0) ^ 0.6 GeV. 
The critical density from pr2T|) . pL22|) is 



p^(criO - L _ }i ) 0.12 GeV^ = 15 fm'^ (11.23) 



so that baryon number density is 



2 



PBar = y ^ 5 /m-3 (11.24) 



This is 30 times normal nuclear density, Uq — ^ fm ^ and the pressure is very high, P^rjj = 125P„„c/(2no). 
It is clear that nuclear matter cannot exist as usual baryon matter for density of 5 baryons in fm^. 
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III. DENSITY EFFECT ON CONFINEMENT POTENTIAL BELOW PHASE TRANSITION 

Consider a white system of a quark and a heavy antiquark Q, which in baryon can be replaced by the string 
junction. Then in the partition function Z = j DADip exTp{ip{D + mq)ip) one can average ever DA and keep the 
quadratic in term in the exponent, which yields Z = J Dipexp J ^{x)ipix)ip{y)ip{y)J{x,y)dxdy.. Here J{x,y) is 
proportional to the quadratic correlator D^{x), J{x, y) = du dvD^{u — v). 

In case of confinement {D^ ^ 0), J{x,x) ^ a; at large x and as shown in [79| this property yields simultaneously 
confinement and chiral symmetry breaking (CSB). As was shown in [t^, the effective quark mass M° (growing at 
large x due to confinement) and quark propagator S are found from the system of equations (we neglect emission of 
pious) 

Mi°)(a;, y) = J^^(x, y)TrSgix, y) (III.25) 



iS,{x,y) - (x|(a + m, + Afio))-i|y). (III.26) 



(0) 



Solutions of pil.251) . (IIII.26P are different for light or heavy quark q. For heavy quark, mq — oo, one neglects M. 
in Sq and has 

Sq{x, y){mq ^ cx)) ~ ^(3) _ y) (jjl.27) 

M^f\x,y) ^ J(x,x) ^a|x-x(0)| (III.28) 

For light quark q one solves the system using relativistic WKB for Dirac particles and again obtains CSB and 
confinement for large x. That was done in [t^ for the case of zero density (vanishing p,q). 

Now we turn to the case of nonzero /iq = /x. As will be seen, for heavy quarks nothing happens -linear confinement 
at all r, unless /i ^ niq (for standard nuclear density /i = 0.3 GeV.). 

Let us study now the case of light quarks, following the paper [80|. The effective light quark mass can be written 

as 

M(x, y; ^) = J(x, y)74A(x, y; ^) (III.29) 

where 74A(x, y, /x) is the time-independent quark Green's function, cf. (jIII.25p . It can be written in terms of one- 
particle eigenfunctions ipn{x) as 



A(x, y; ^) = ^ Vn(x)sign(e„ - Ai)-0+(y) 



= Ao(x,y)- AA(x,y) (III.30) 

In [131 V'n(^) found in the relativistic WKB method and the approximate summation over n can be performed, 
yielding for Ao(x, y) the smeared (5-function form for large |x|, |y| and |x — y| ^ |x|,y|. In [s^] also the term 
AA(x, y) in (jlll.BOp was calculated in the same way, which is nonzero for /i 7^ 0, and distorts the linear confinement 
at distances r < ^/cr, where a = 0.18 GeV^ is the standard string tension. The resulting picture for the averaged 
M(r) = J (i^r'M(r, r', /x) is given in Fig|31 where M(r) is split into a Lorentz scalar and vector parts, 

Mscai = (Tr9{ar - ^i), NUect{r) = -2(jre{^i - ar).) (III.31) 

In MQS we shall for simplicity consider all quarks attached to the same string junction (this is not crucial for quarks 
in the modified confining potential, since near the string junction the scalar part vanishes, but the vector part can be 
rewritten in the form of pairwise potentials [8l|). Consider now a system of N quarks around one string junction in 
the nuclear medium of ordinary nucleons. Using the local limit for interactions of a light quark with string junction, 
Eq. (|III.3ip . one can calculate the change of masses of white states due to density. 
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FIG. 3: Interactions Mgcaiir) (solid line) and Mvect{r) (dashed line) as functions of distance r, with b = /^/o". Dotted line and 
dashed line show the qualitative smoothed form of both terms respectively. 




Assuming that multiquark states located at different string junction points do not overlap (see below the note on 
states radii), one can associate the parameter /i for a given quark with the largest value of £„, occupied by other 
quarks, belonging to the same string junction point. In what follows we shall show, that in some situations it will 
be advantageous for two or more nucleons (with different string junctions) to coalesce into a common state with one 
string junction, building in this way the contracted potential (|III.31I) . 

Note, that for one nucleon (3 quarks) the contraction mechanism does not work and fi can be taken at zero value, 
since no forbidden states exist for each quark of a given color. This, again, is true, provided one-quark states do not 
overlap at the given density. The situation changes, however, for two nucleons, since e.g. one can add to the first 
nucleon, i.e. 3 quarks in S-states, 3 quarks in P-states relative to the same string junction point. The latter will move 
in the scalar part of the potential A(x,y), with ^ being equal to the S-state Dirac eigenvalue £5, thus the energy 
ep{p- = es) may become less than £5(^1 = 0): 

ep{^l = es) <es{ii^Q) (III.32) 

The same type of inequality may occur for higher L states, when more nucleons coalesce to the same Multiquark 
State. This would imply instability of nuclei with respect to transition to MQS matter. 

Note, that /is 7^ 3/i, since /i is attributed to MQS string junction and does not grow with unless different MQS 
start to overlap. Hence for /xs > (jIII.32|) should be replaced by 

ep(M = £s) < es(M = 0) + £B(Fermi). (III.33) 

To check pil.33p in [sOl £s and ep were computed numerically with the potential (|III.31|) . and in perturbation theory, 
considering the change of potential in (jIII.31[) as a small perturbation. The results are given in Tables IIIIIIII The 
calculated quantities are AM5 = 3A£„(5) and AMp{3q) = 2A£„(5) + A£(P), where A£„ = £ri(A*) ^ £ri(0). One 
should note at this point, that the vector part of interaction written in (IIII.3ip for the case of light quark with heavy 
antiquark, and for the case of 3 quarks interacting with string junction the vector part is transformed as an addition 
to the color Coulomb potential with coefficient 1/2 (see [8l| for explicit derivation of scalar and vector interaction in 
the nucleon). This prescription was used for our calculations. 

From these numbers one can obtain es{fJ.q = 0) « 0.52 GeV, ep{^q = 0) = 0.77 GeV, ep{^q ~ 0.52 GeV) = 0.57 
GeV. This is close to the es{nq = 0), i.e. one can expect the formation of MQS for eb (Fermi)=50 MeV, which 
corresponds to critical density p = 3pQ. 

Also, the average radii of S- and P-wave quark states were computed (see Tables Hilllll ) and it is found the radius is 
decreasing with a rise of fi, at least for the moderate p values, so it is reasonable to neglect in the first approximation 
the possible overlapping of multiquark state wavefunctions attached to different string junction points. 

As we shall discuss below, this might be important for new heavy- ion colliders (FAIR, NIC A,...), and crucial for 
neutron stars and supernovae. However, we have still neglected interaction between baryons in nuclear matter. In the 
next sections we show that this interaction can be described in terms of MQS formation in the intermediate states, 
and a new theory of nuclear matter based on MQS can be formulated. 
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TABLE II: Mass shifts of 3q system with aU quarks in the S-state (purely numericaUy, and perturbatively, for comparison) and 
the S-wave average radius, due to nonzero fiq 



/i, MeV 


AMexact, MeV 


AMpert., MeV 


(r)s, fm 


100 


-5.8 


-5.7 


0.505 


200 


-75 


-70 


0.487 


300 


-290 


-265 


0.460 


400 


-630 


-603 


0.464 


500 


-1008 


-1020 


0.501 



TABLE III: Mass shifts of 3q system with two quarks in the S'-state and one quark in the P-state (purely numerically, and 
perturbatively, for comparison), and the P-wave average radius, due to nonzero /iq 



M, MeV 


AMoxact, MeV 


AMpert., MeV 


(r)p, fm 


100 


-4.0 


-3.9 


0.613 


200 


-53 


-50 


0.605 


300 


-221 


-199 


0.558 


400 


-530 


-481 


0.473 


500 


-917 


-873 


0.448 



IV. QUARK COMPOUND BAG MODEL 



In this section we describe the general two-channel formalisni of the nucleon-nucleon interaction that provides 
dynamical framework for the P-matrix formalism and is used in Sect. VI to build the relativistic formalism and to 
construct in Sect. X the Bethe-Goldstone G operator for nucleons in nuclear matter. 

For simplicity we consider a hadronic channel consisting of two spinless nonrelativistic hadrons, which is coupled 
to another channel, where any number of quarks and gluons are present. We denote the hadronic channel with a 
subscript h and the quark channel with a subscript q. The wave function is a column with two components (^E'/j, ^i^), 
while the interaction term in the Hamiltonian is a 2 x 2 matrix. 

^'^'^ ^ . (IV.l) 

We do not specify the quark dynamics and even the type of variables, on which Vqq is acting. At this point we can take 
the description of the quark- gluon system as general as possible and we shall use only the property of confinement, 
i.e. that the eigenstates of the quark Hamiltonian are discrete states confined to a finite region in space whatever the 
binding energy is. 

The coupled dynamical equations for "^h, ^ q are 

{TH + Vhh-E)-^h = -Vhq^q, (IV.2) 

Lq^q + Vqh'^h = (IV.3) 

where Lq is some quark-gluon operator (its exact form is irrelevant for our purposes) and is the kinetic operator 
for hadrons and E is the full energy, E = y/s; Vhh is in general a nonlocal energy-dependent operator. Th in general 
is a relativistic operator. 

One can solve formally equation (jIV.3[) to express in terms oi 'i'h- 

^q = -GqVqh'fh (IV.4) 

where the quark Green function has a spectral representation 
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The quark eigenstates ^'^ and eigenvalues E^, belong only to a discrete spectrum because of the confinement. Using 
(|IV.4p one can exclude the quark channel and rewrite the hadronic equation pV.2p as 

{Th + Vhh + Vhqh)^h = (IV.6) 

where 



Let us discuss the structure of the quark-induced hadronic interaction Vtqh- 

1. It does not contain any quark or gluon degrees of freedom, which are present in and Vhq- Indeed, in (|IV.7I) the 
summation over all quark degrees of freedom is implied by the repeated subscript g, so that the only remaining 
coordinates are hadronic relative distances, entering via Vhq and Vqh- 

2. Vhqh is separable Hermitian interaction; if we denote 

U{r) = Vk,^\. (IV.8) 

Vhqh has the form 

Vhqh^l^ • (IV.9) 

V 

3. The numerator of pV.9p depends in general case on energy E through the quark-hadron interaction Vqh- The 
most characteristic energy dependence enters through the denominator in (jIV.9|) : for E near E^, the quark- 
induced potential Vhqh can be infinitely large. 

4. Since quarks are not allowed outside the bag the function should be nonzero only inside the bag. Hence 
it can be easily approximated there by any series on a finite interval of r. 

We assume now that the solution of the hadronic problem (jIV.6p without Vhqh is known; the corresponding regular 
at origin scattering wave functions and the Green function we denote by X{r) and Gh respectively. In terms of these 
pure hadronic quantities we can rewrite the full coupled quark-hadron equation (|IV.6P as 

^h{r) ^ X{r) + ^ (GhMif.^h) ^j^^Q^ 

1/ ^ 

(Ghfu) = j Gh{r,r')Ur')dr' 

U.^h) = / ^h{r)f;{r)dr. (IV.ll) 



where we have used the notation 



The constants {f^^h) are to be found from the system of linear algebraic equations 



where 

iUGhf^)^ I I f:{r)Gh{r,r')f,{r')drdr'. 



The resulting picture is particularly simple when E is close to some eigenvalue E^^ and all other terms in the sum 
in (jIV.lOl) with v ^ [l can be neglected. In this case we have 
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At this point it is possible to discuss the general properties of our solution ^'/i(r). First of all, multiplying both 
sides of pV.13p with f*{r) and integrating over dr, we obtain 

iU^h)^0, for E = Ef,. (IV.14) 

This property does not mean that vanishes inside the bag, where /^(r) is nonzero. Indeed, ^'/i(r), being complex, 
may have no zeros at all inside the bag, still the condition (jIV.14|) can be satisfied. One example of this kind obtains 
if we take /^(r) — 9{R — r) and Vhh = 0. No definite consequences can be deduced for the observables in that case. 
At this point we make an ansatz, that hadrons are coupled to the QCB degrees of freedom only on the surface of 
the compound bag. As we shall see below, this assumption immediately leads to explicit expressions for the wave 
function and the observables. It is crucial for all results of the paper. Moreover we show below that this assumption 
leads naturally to the P matrix formalism of Jaffe and Low [50| and therefore experimental evidence in favor of that 
formalism presented in [50] justifies our assumption. Moreover in Refs. (49l . [5l| additional experimental evidence is 
presented in favor of the stated ansatz. 

From theoretical point of view the assumption, that the hadrons effectively are not coupled with the bag constituents 
inside the bag, but only at the surface, is not new. In the framework of the nuclear cluster models, (see the book (82j ) 
it is well known that the region near the surface of the compound nucleus plays the crucial role. 

In [s^ an additional support is given to the assumption in terms of quark cluster models. The general character 
of these considerations allows to apply it also to the quark systems. Accordingly we assume the following simplified 
form for /^(r): 

}^{r)=c^5{r~b) (IV.15) 

where h is connected to the QCB radius R and is some constant. With the choice (jIV.15[) the solution '^h{f) 
simplifies: 

^.ir)^Xir)^^^4Sm_ (IV.16) 

where 7^ = |c^ p. 

From (jIV.16[) follows an important property of '^h{r), namely 

= 0, for r<h and E ^ E,,, (IV.17) 

since 

G;,(r,/)=X(r<)y(r>) (IV.18) 

where Y{r) is a nonvanishing at r = solution of the purely hadronic equation (|IV.2p with Vhq — 0. It is convenient 
to choose X{r),Y{r) with the asymptotics X{r)r^oo ~ sin(fcr+jo(fc)) 

Y{r)r^oo - 2/iexp[i(fcr + ,5o(fc))], (IV.19) 

where /i is the reduced mass (reduced energy in relativistic case). 
Let us define the factor F{E) as follows 

|*„(r)|2 =F(£;)|X(r)|2, r < b. (IV.20) 

F{E) shows how much 5'/i(r) decreases inside the bag because of the presence of the quark states; exact form of F{E) 
is: 



F{E)^ 



E^-E 



E^ ~ E ~ -f^Ghib, b) 



(IV.21) 



From (jIV.16[) we immediately conclude that the logarithmic derivative of '^h{f) at r = 6 should be meromorphic in 
E, namely 

+ (IV22) 
To obtain (IIV.22|) we use the property of the solution X, Y: X{r)Y'{r) - X'{r)Y{r) = -2^. 
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An important quantity entering (IIV.2ip and (|1V.16|) . Gh{b, b), can be expressed through Pq = ^(bj ^^^^ latter can 
be extracted from the experimental data, see [Iq!]). Assuming that there is no interaction beyond r = b, we can use 
PV.19|) as the exact values of X{r) and Y{r) for r > b, and obtain from (|IV.18P : 

Gh(b, 6) = ^ sin^ikb + So)(cot(kb + 5q) + i). (IV.23) 
k 

On the other hand we have Pq = kcot{kb + Sq) so that 

Gkib,b)^-^. (IV.24) 

Pf) — ik 

It is instructive to evaluate ''y^Gh{b,b) which enters (jIV.16|) . For NN scattering from the analysis of [1^ we have 
Po « 0.4 GeV, 7^ « 0.1 GeV. Therefore e.g. at £; = 0.1 GeV, -f^Gh{b, b) = (0.15 + i0.12) GeV. 

The ratio F{E)/F{0) is effectively zero for positive energy, e.g. for E — 0.05 GeV this ratio is ^ 0.01. Let us now 
consider many quark eigenstates instead of one. 

In this case the general result (IIV.10[) for is valid where {f^'^h) are to be found from the solution of the system 
of equations (|IV.12p . The results simplify very much if we use the ansatz analogous to (|IV.15P : 

U{r)^cJ{r-b), v^l,2,...N (IV.25) 

with the same range b for all quark eigenstates. 

The ansatz (|IV.25p can be approximately true for those bag eigenstates where the number of quarks and energy 
eigenvalues do not change much. 

Taking into account (IIV.25|) we obtain {f^giiftj,) — c*Gh{b,b) and [f^'^h) — c*^/i(6). Inserting these values we 
rewrite Eq. (jlV.lOp as follows: 

Gh{r, b)X(b)y^,^-^ 
= X{r) + / : X (IV.26) 

1-Gh{b, b)J2^E^ 

where -y^, = \c^ \'^. 

From pV.26|) one can see that whenever E approaches some eigenvalue E^, the hadronic wave function vanishes for 
r < &. In this way the property (jIV.171) generalizes to the many quark eigenstates case: the form of the wave function 
(IIV.26I) is similar to the one-state case (jIV.lGp if in the latter one replaces in the denominator by the quantity 

X -1 

E^-E ) • Under the same condition as in (|IV.25P we obtain the following expression for P matrix 

P = k cot(fc5 + (5o(fc)) + Y ^^^^ (IV.27) 
^ E-E, 

and all arguments about the correspondence of the quark eigenstates ( "the primitives" the terminology of Jaffe and 
Low [50| ) and poles of the P matrix remain true. Here (5o(fc) is the phase due to interaction Vhh which is assumed to 
vanish for r > b. 

In practical calculations one is using the potential of the form [i^, [s^ 

V - Vhqh + VhhO{r - b) (IV.28) 
and Vhqh for arbitrary NN angular momentum I is 



T/ (r-r.^^ ^ \-V* mr,E)Ur',E) ^^ 

Vhqh{Y, I- ) = ^ 2^ Yj^i^ ^ Y.^is, (IV.29) 



n{T,E)Ur',E) 

m 

while fiy{r, E) can be represented as 



U{r,E)^-cJ{r-b)+x,{E,^E)7^,{r), (IV.30) 
where r]i,{r) describe the relative NN motion inside the MQS' and are chosen in the form 



Tlu{r) = iV.Jz(/3.r)/3,r, / r,t{r)dr = 1. (IV.31) 
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TABLE IV: Even and odd operators O. 



Even 


Odd 


even 


odd 


even-odd 


odd-odd-odd 


1 

ns 


7m 


r"i75 


a ■ ■"W' 






{ d + d )i,(7^i^ 




(d + d). 


i d fj. 


(9 + 9 )^i750-,i„ 





Finally, all coefficients are to be fitted to the experimental phases, usually in the interval [ 0; 500 800 MeV ]. As 
a first approximation we shall assume, that Vhh = and only Vhqh is retained; as is shown in [i^, [5l| in this way 
one obtains a good fit to experimental NN S phases for energies up to ^ 500 MeV. In Fig. [3] comparison to the 
experimental phases is given for the ^Sq and '^Si states obtained from relativistic-expressions of Sect. VII. 

We thus introduced in the NN system the interaction due to MQS PV.29I) . which can alone (or in combination 
with usual OBEP) describe NN force, and the same in 3A^, AN... systems. At this stage /„,£"„ are found phe- 
nomenologically from NN phases, and due to |40|] are changing with matter density. One can use Vhqs changing with 
matter density. One can use Vhqs both to reactions. At the next stage one should calculate fm En microscopically 
from QCD Lagrangian (e.g. in the framework of the FCM). 



V. EFFECTIVE LAGRANGIAN OF QUARK COMPOUND BAG MODEL 

In this section, we construct effective relativistic QCB Lagrangian for the description of nucleon-nucleon interaction. 
Two-nucleon currents entering the effective Lagrangian are constructed as bilinear combinations of the nucleon 
wave functions: 

J{x) = vp^(a;)0^'(a;). (V.l) 

The nucleon wave functions, ^'(a;), carry bispinor and isospin indices, 5'c(a;) — C'^'^{x) is C-conjugated wave function, 
= (*c) = — ^'"^C, and C = ir^ij'^j^ in the standard representation [83|. The currents (jV.ip annihilate dibaryons. 
The creation currents have the form 

J+{x) = §(a;)0*c(a;). (V.2) 

The matrix O admits expansion over the Dirac 7-matrices and the Pauli r"-matrices. Many terms of the formal 
expansion vanish, since 4'(x) anticommute. Under the permutation of two nucleon fields J{x) transforms identically 
provided 

C'^O^C = -O. (V.3) 

Transformation properties of the operators O and the nucleon fields under the C-conjugation are given in Appendix 
B. 

Matrices O entering the effective Lagrangian are therefore antisymmetric (odd) under the C-conjugation. In Table 
IIVI we show elementary even and odd matrices and the composite odd ones which obey (IV. 31) . They are combined 

further with the first-order differential operators i d ^ = i{— d + "^)^ (odd one) and {d + 9 )p (even one). 

The matrices t"757^ and T^a^^i d y are not shown, since currents they produce with free nucleons represent full 
derivatives: 

^'ct"757a.* = -a^(*cT"i75*), (V.4) 
ra 

^cT^a^yil^y^ = 2a^(^,T"*). (V.5) 

Full derivatives shifted to the dibaryon vector fields give vanishing contributions to the Lagrangian. 

We restrict ourselves with scalar and vector dibaryon currents of positive and negative parity. This is sufficient to 
provide the phenomenological description of the NN scattering in J = 0^, 1^ channels {S- and P-wave scattering). 

The list of possible two-nucleon currents is shown in Table |Vl 
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TABLE V: Two-nucleon currents, associated dibaryon fields, and their quantum numbers. 



V_/ Ul I UllL 


1 11 Q "Vi rj^TI 

J_/iL)cLI y Uli 


\^^-J ) 






(1,0-) 








*cj d ^* 


X+M 


(0,1+) 




X+M 


(0,1+) 


X-M 


(0,1-) 


a. (denser""*) 

*cr"j75Cr'"'iV„* 


X-M 


(0,1-) 


X— /.^ 


(1,1-) 



A third vector current 4'c7^5' with quantum numbers (0, 1+) can be expressed in terms of the two ones hsted in 
Table fVl using the Gordon's expansion 



1 - ^ 
2m 



2m 



The vector dibaryons are characterized by two coupling constants. Recall that photon couplings with the nucleons 
are characterized by two coupling constants also. 

In the channel (0, 1+), the pseudovector coupling with the current d ^4* gives vanishing contribution at two- 
nucleon threshold, so the 5- wave NN scattering length is determined completely by the pseudotensor coupling (0,1+). 
In what follows, we neglect by the pseudovector coupling. 

In the channel (0, 1-), a linear combination of two currents enters the P-wave NN scattering length, so we set the 
vector coupling equal to zero and redefine the tensor coupling. 

Finally, every dibaryon is coupled to one two-nucleon current. 

The dibaryon fields Lpf,, XP^t ^-nd x'p^ carry (suppressed) indices of radial excitations, index P = ±1 stands for the 
parity, a is isospin index, and /i is the Lorentz index: 



The coupling constants are defined accordingly: 





\ 






\ 
















, XiJ. 
















V 


) 









\-fJ- 
X2+I1 



(91-] 




( hi+ 






(K+ 


91+ 




hi. 








92- 




h2+ 






h'2+ 


92+ 




h2- 






h'2- 


V- ) 




['■■ 


J 




['■■ 



(V.6) 



The effective Lagrangian splits into free and interaction parts: 

r - rlo] + rW , /-M 

— '-int ' '-inf 



'^Int describes interaction of dibaryons with nucleons, while C^^}^^ describes a four-fermion contact interaction. 



The free part has the form 



+ a^^"+9^(^" - (^"+m2o< 



£[0] = *(iV-m)4' 

duX'j.dvXii + X^ AfgiXM 
duxTd^xZ+xT^liXl 



(V.7) 



(V.8) 
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The interaction parts are as follows 



'[1] _ 



[2] 



h/_ h/+_ 
2m 2m 
= Aio(*r"r*,)(*,r"r*) 
+ Aoi(*ra^.*e)(*cra^.*) 



(V.9) 



(V.IO) 



Here, m is the nucleon mass, Afjj are dibaryon masses (matrices) diagonal in the radial and parity quantum numbers 
and 



/ 1 



10 



v; ; 



The coupling constants g, /i, /i', and A/j depend on the radial and parity quantum numbers as indicated e.g. by 
Eq. (jV.6p . C is Hermitian by the construction. 

There exists six various vertices for the Lorentz structures: {1, 7^t, o'/ji^j ^757 ^Tsc^ti/}- The vector and tensor 
ones correspond to two spin-1 fields, vector and pseudovector ones. So, we have at most two spin-0 fields with parity 
P = ±1 and two spin-1 fields with parity P = ±1, altogether four fields. The isospin degrees of freedom provide two 
additional structures, so we might have eight dibaryons. However, not all of them are coupled to the NN channel. 
The system of two nucleons may carry quantum numbers shown in Table IVII 



TABLE VI: Two-nucleon states with lowest quantum numbers. 









(0,1+) 


^So 


(1,0+) 


^Pi 


(1,1-) 


^Pi 


(0,1-) 


^Po 


(1,0-) 



We thus have five possible dibaryons (/, J^) = (1, 0^), (0, 1^), (1, 1^). Missing are three dibaryons with the exotic 
quantum numbers (0, 0^), (1, 1+), that do not correspond to two-nucleon states. From other hand, states (0, 0^) have 
the vanishing couplings with the NN channel, since = ^c^l^"^ = 0. In the (1, 1+) channel one has, in addition, 
^ct"7m* = *cr"(T^.* = 0. 

The (1, 1^) channel is not exotic. It is included into the effective Lagrangian through the vertex ^cT^i^^cr^i, d . 
The vertex ^cT'^a^^ d is full derivative according to Eq. (jV.Sp . 



VI. IN-MEDIUM NUCLEON AND DIBARYON PROPAGATORS 



We consider symmetric nuclear matter with equal proton and neutron fractions. 



A. Nucleon propagators in ideal Fermi gas 



The in-medium nucleon propagator is defined by 

iSpix) =< T«'(x)*(0) > 



(VI.l) 
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The propagator depends on the Fermi momentum pp . In the momentum representation, 

Spij)) = j dxe'P'^SFix) (VI.2) 

We need the in-medium nucleon propagator constructed out of the C-conjugated nuclcon fields 

iS'pix) r*c(a;)*c(0) > . (VI.3) 

The momentum representation of Sp{x) is derived from Eq. (|VI.2|) . The propagators are related by 

S-pix) - C^Sl{-x)C, (VI.4) 

S'p{p) = C^Sl{^p)C. (VI.5) 

The C-conjugation matrix is defined in Appendix A. 

In the momentum representation, the propagators of the ideal gas have the form 

Sf{p) = SUp) = (VI.6) 
p — m 

The signs of imaginary parts of the pole positions on the complex energy plane should be specified separately. 



The poles of Sf(p) approach the real axis from the upper half of the complex po-plane for SRpo < = +\/"?+p|^ 
and from the lower half of the complex po-plane for 5Rpo > Sp{p) has poles shifted to the upper half-plane for 
Sftpo < ~M ^-i^d to the lower half-plane otherwise. 

Sf{x) describes therefore propagation of nucleons in presence of the Fermi sphere made up of nucleons. Sp{x) 
describes propagation of nucleons in presence of the Fermi sphere, however, made up of antinucleons, the Fermi 
momentum being the same. 

The plane waves expansion of the in-medium propagators has the form 



iSf{x) = I [e-^^'^A+lp) {e{xo)e{\p\ - pf) - 9{-xo)e{pF - IpD) + e'f^A_(p)0(-xo)] , (VI.7) 



^S'f(.^) = I [e^P^A^ip] iei^xo)9i\p\ ^pf) - e{x^)e{pF - \p\)) + e-^P-K+{p)e{x^)\ , (VI.8) 



where p — {E{p), p), E{p) — +1/771^ + p-^, and 



are projection operators that obey 



C^A±(p)^C = A^(p). (VI.IO) 

B. Dibaryon propagators 

The dibaryon propagators with quantum numbers (1, 0^), (0, 1^), (1, 1^) are defined, respectively, by 

iA'°'^{x~y) = <Tip"{x)ip^+{y) >, (VI.ll) 

iD'^^ix-y) ^ <Tx^.{x)xt{y)>, (VI.12) 

iD'^fix - y) = < rx^(x)xe+(y) > . (VI.13) 

The radial indices are suppressed. The coupling to the NN channel produces mixing of the radial excitations. 
The free propagators are given by 

^"^(P) = ^^^^"'^ (VI.14) 
P 10 

DM - (VI.15) 

D^^ip) - (V1.16) 

p — Mil 
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Tensor structures of the dressed propagators can be factored out 

A"'P{p) = S^^A'ip), (VI. 17) 

D^P) = i-9t..+P^.P./p')D'{p), (VI.18) 

D'^P{p) = 6^f'{-g^.+p^pJp')D[{p). (VI.19) 

C. In- medium dispersion law for interacting nucleons 

In systems with interaction, self-energy operators of fermions appear. They can be expanded over the 7-niatrices: 

t[p) = tv{p) + ^s{p). (VI.20) 
t%p) = + (VI.21) 

where = Tjy{p)^^ and = (p)7^. The C-conjugation gives (cf. (|VI.5|) ) 

= C^t^{-p)C. (VI.22) 

We thus obtain 

S^^(p) = (VI.23) 
- Ss(-p), (VI.24) 

and 

^^''(p) = 75S(-p)75- 

The diagram representation of the self-energy operator arising due to the dibaryon-induced nucleon-nucleon inter- 
action is shown in Fig. Its non-relativistic version is C/2(p,p) defined by Eq. (IIX.17p . Equation (jlX.lSp is the 
analogue of the one-loop approximation shown in Fig. |4l 




S(P,Pf) = 

FIG. 4: Nucleoli self-energy operator. The double solid line with the shaded blob shows the dressed dibaryon propagator, the 
single solid line the shaded blob shows the dressed nucleon propagator. 

The in-medium propagators have the form 

1 p + TO + p) 



S'f{p) = 
S'^ip) = 



p-m-t(p) Hp) 

1 _ p + m + S(— p) 

p~m~t''{p) <l>c{p) 



where 

= {p-^v{p)r-{m+^s{p)r, 

Up) = {p-T.<y{p)f-{ra + i:%{p))\ 
Using Eqs. (|VI.23|) and (|VI.24|) we find 

Up) = c^{-p). (VI.25) 

The dispersion laws of elementary excitations of nucleons are determined from equation 

(j){p) = 0. (VI.26) 

This equation gives two solutions po = i£'±(p) with positive and negative energies of the excitations. Respectively, 
positive and negative energy solutions of equation cj)c{p) — are po — ±£'=f(p)- the ideal Fermi gas, i?±(p) — 

+ \/'W? + p2. 

The plane waves expansion of the in-medium propagators has the form 
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^S'A^) = y'^[e-^^'+-AV(p)(0(a;oM|p|-pi.)-0(-xo)%i.-|p|)) + e^f--A'_(p)(?(-xo)], (VI.27) 



dp 

where p± — {E± (p) , p) and 



'Spci^) - / (2^)3 [e''+"AV^(P) {Oi~xo)e{\p\-pF) - eixo)e{pF - |p|)) + e-^P-^A'l{p)6{xo)] , (VI.28) 



with = d(f){p)/dpQ and = d(j)c{p)/dpo. In the ideal Fermi gas, Eqs. (IVL7P and (jVLSp are recovered. 



D. Dyson equation for nucleon propagator 

is calculated to one loop, contributions of the antinucleons are neglected. The loops are formed thereby by 
nucleons and nucleon holes in the Fermi sphere. Such approximation is close to the mean-field approximation. 
Shown in Fig. [5] is the system of Dyson equations, 

S'p = Sf + Sf^S'p, (VI.31) 

for the in-medium nucleon propagator. The virtual dibaryons in the loop contribute to the self-energy operator. The 
four-fermion interaction will be included later. 
Starting from Lagrangian (jV.7|) one gets 



\±ip) = I ^g+r'^^A'"^ip-p,)^SMr^^g 



+ 



(2^) 

d^Pc /i+ h 

d^Vr h'+ 



h' 



in \4 r, i<^p.T{P + Pc)TT'^ilzDl2^{p~Pc)iSp{pc)TPi"tzi(J^,„{-pc-~p)a7^- 



(VI.32) 



In Eq. (|VI.32p summation over the radial numbers and parity states (F = 1,175) is assumed. The integrals over the 
timelike component of Pc can be calculated assuming (f>{p) has simple zeros in the complex energy plane at po = 
±£'-l-(p). The negative energy states are neglected. The remaining three-dimensional integrals describe contributions 
to the self-energy operators from the nucleon holes. 

The self-energy can be expanded over the different channels 



±ip) = j2^jjHp)- 



(VI.33) 



1.7'= 



The separate contributions for (/, J^) = (1, 0=^), (0, 1±), (1, 1") are given in Appendix B. 



o 



o + 0- 




o 



FIG. 5: System of Dyson equations for nucleon propagator S'p{p). The single and double solid lines show the nucleon and 
dibaryon propagators. The shaded blobs indicate the dressed propagators. Summation over dibaryons in the loop is assumed. 
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E. Bethe-Brueckner equation for dibaryon propagator 



Dibaryon self-energy operator to one loop is determined by two-body nucleon Green function, its diagram represen- 
tation is shown in Fig. [51 The system of equations of Fig. [5] is equivalent to Eqs. (lIX.ip which determine two-nucleon 
scattering in the vacuum (T-matrix) and in the medium (G-matrix). The potential V2 is the analogue of the dibaryon 
propagator in the vacuum, G2{W) is the analogue of the dressed dibaryon propagator, (^262^ is the dibaryon self 
energy shown in Fig. [S) 

The corresponding analytical expressions have the form 

/j4 
-^gSp[T^rs'Ap + Pc,PF)r''rs'p,ipc,PF)]9+, 

^yt^iP.Pp) = -2i / -77-^ — S^[{-^(T^aPa)yS'p{p^Pc:PF)i(T^,TPT^S'p^{pc,PF)]-z—, 
J (ztt)* 2m 2m 

^^;^^Sp[(-T^i75ZCTj,„(p + 2pc)CT)5'^b + Pc,PF)«(^,<r(p + 2pJ^T"i755^c(Pc,PF)]^- 

Here, Sp[...] is taken over isospin and bispinor indices. The self-energy operators are diagonal in (/, J^) and represent 
matrices in the radial numbers. 



n(p,pp) = 



FIG. 6: Dibaryon self-energy operator. The single solid lines with the shaded blobs are the dressed nucleon propagators. The 
loop produces mixing between dibaryons with different radial quantum numbers. 



The system of Dyson equations for the dibaryon propagators has the form: 

A'"'^(p) = A"^(p) + A""(p)n'^^(p)A"'^(p), 

D'^M - D^,{p) + D^r{p)nra{p)DUp), 

KfiP) = DfAp)+D"^Up)Uil{p)D'Z^ip). 
These equations are depicted in Fig. [71 




F. Link to optical potential model 

Probing particles acquire new decay channels in the medium, connected to collisions with the surrounding particles, 
and become unstable. Phenomenologically, the in-medium spin-1/2 particles are described by the Dirac equation 

(iV - mi + ^F)* = 0, (VL34) 

where F is collision width 

F = pavj. (VF35) 

Here, v is velocity and 7 Lorentz factor. The cross section can be expressed in terms of the forward scattering 
amplitude with the help of the optical theorem 

a=^lmf{0^O). (VF36) 

p* 
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The self-energy operator and the forward scattering ampUtude are analytical functions. Given the imaginary part 
is known, the real part can be restored using the dispersion relations. If, however, the scattering amplitude is known, 
one can write 

m = -P^^fiO - 0). (VI.37) 

TOir7l2 

where m2 is mass of the particles of the medium. 

Near the threshold, one can expand the amplitude over the partial waves and keep the lowest order L = terms: 
f{0 = 0) w e'^'^'^ smS{s)/p* where S{s) is the S -wave phase shift. If probing particle is identical with particles of the 
medium, one has nii — TO2 — m, f{6 = 0) « 2e^^^'^'> smS{s)/p* and 

= _p-^e**(^) sm5{s). (VI.38) 

The symmetric nuclear matter consists of the two types of the particles, so the statistical factors modify Eq. (IVI.38I) . 

Let a spin-up proton scatters on the in-medium protons in the (1,0"*") channel. The in-medium protons have 
spin down. The statistical weight of the configuration in the S = spin wave function is 1/2. The proton-proton 
collisions contribute to the self-energy operator with the weight 1/2 x p/A, where p/i is the spin-down proton density. A 
spin- up proton scatters on the in- medium neutrons also, with the weight factor 1/2 for the 1=1 channel and the weight 
factor 1/2 for the spin-down neutrons. The spin-isospin factor in the 1 = 1 channel becomes 1/4 x 1/2 x 1/2 = 1/16. 
The in-medium protons and neutrons contribute altogether with the factor 1/8 + 1/16 = 3/16. We therefore replace 
p 3p/16 in Eq. (|VI.38|) . The (0, 1+) channel can be analyzed in the similar fashion. 

From other hand, in the low-density limit of the symmetric matter, Eqs. (jVI.32p give 

S/jp(rt = -p-^e'^./^(^) sin S J jp{s)^^^^ + ..., 



where (/, J^) = (1, 0+) and (0, 1+) in agreement with the optical potential model. 

It is instructive to compare the nucleon self-energy with the Relativistic Mean Field (RMF) model [l^,!!^, according 
to which 

The first term describes the w-meson exchange between the nucleons, the second one originates from the cr-meson 
exchange. In the RMF models ^ 7" — 1 (numerically). The strong cancellation of the large vector and the large 
scalar mean fields is the typical feature of the RMF models. 

In our model, the nucleon self-energy in the low density limit is proportional to 70 + 1. The first term can be 
interpreted as the vector mean field, the second one as the scalar mean field. The mean fields are of the same sign 
and of the same order. 



VII. DIBARYONS IN PHASE SHIFTS ANALYSIS OF NUCLEON-NUCLEON SCATTERING 



In order to describe the in-medium properties of nucleons and dibaryons, one has to fix parameters entering the 
effective Lagrangian. Those parameters determine the NN scattering phase shifts which are known experimentally. 



A. Dibaryon self-energy operator in the vacuum 



The dibaryon self-energy operators are matrices proportional gg^ , hh^ , and h'h'^ . Such a circumstance allows to 
factorize the radial numbers and the parity structures. Also, the isospin and Lorentz structures of the self-energy 
operators can be factored out: 

n^-{p) = s^"gn,opip)g+, (vii.i) 
Ti.M = {-g.^ + ^)hnoiP{p)h+, (VII.2) 

nf^(p) = 6^"{-g,^ + ^)h'U,,-{p)h'+, (VII.3) 



where 
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vij^iqcPc) 



i2n) 

The functions ipjjp{qc,Pc) have the form 
2i 



{ql — m? + — m? + iO) 



(VII.4) 



'4^wp{qc,Pc) 
^oip{qc,Pc) 
<Pii-{qc,Pc) 



Sp[T"r(ge + m)T"r(-p, + m)], 



3 

2i 1 

3 Am? 

2Z 1 PvP L 

-r-. — ^(-Si'M + ^V^)Sp[r"i75(-ia-^cr('7c - Pc)a){qc + m)ia^i,r{qc - Pc)TT"i-f5{-pc + m)]. 
9 4to^ 



(-.9i.M + ^^)Sp[r(-zcrt.crPCT)(gc + m)i(Tf_,rPT'r{-Pc + m)], 




FIG. 7: Systems of Dyson equations for the dibaryon propagators shown by the double solid lines. The single solid lines show 
the nucleon propagator. The shaded blobs indicate dressed propagators. 



The discontinuity over the cut (4m^, +oo) can be found to be 

Ilijp{s + iO)-nijp{s-iO) = -~^<^>2{s)(pijp{qc,Pc), 

where qc and Pc are the on-mass-shell nucleon momenta and $2(5) is two-body phase space 

7rp*(s) 



$2(S) 



where p*{s) = y^s/A — m?. In what follows, p*{s + iQ) = k > Q above the threshold and p*{s — iO) = —k. 

The functions fijp{qc,Pc) are given in Table IVlIl They show the correct threshold behavior ^ k^^ . The imaginary 
part of II/jp (s) is negative and corresponds to positive transition probabilities. 



TABLE VII: ^ijp{qc,Pc) and K/jp is different channels. 









'Si (0,1+) 


8js(l -f p*V(3m2)) 


■Kb/{2m) 


'So (1,0+) 


Sis 


nb/{2m) 


'Pi (1,1-) 


64ip*V3 


nmb^/ 12 


'Pi (0,1-) 


16ip*2s/(3m2) 


TTjnb^l 12 


'Po (1,0-) 




Trmb'/lS 



B. S matrix 

We consider first the scalar dibaryons. The dressed propagator can be written as follows 

A'"'5(p) = 5"^A'(p). 

In the one-loop approximation, the system of Dyson equations can be solved to give 

A(p),gnio-(p)5+A(p) 



A'(p) - A(p) 



l~ K,^p{p)IV^op{p) 
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Here, A(p) is the bare dibaryon propagator diagonal in the radial numbers, 



A-V(s)=g+A(p)5 = g+—}-^g (VIL5) 



1 



5r, 



where s = , the sum is taken over the radial numbers. 

The model described in Sect. V is the relativistic extension of the Quark Compound Bag model [4^]. Aq^p(s) is, in 
particular, the function defined in Eq. (5) of Ref. [i^. The poles of A.~^p{s) correspond to the compound Qq states 
i.e. dibaryons, whereas poles of Aqip(s) are the Castillejo-Dalitz^Dyson (CDD) poles [1^1 • The model described in 
Sect. V can also be considered as a version of the Dyson model |865, which is, in turn, a version of the Lee model [87| . 
The four-point Green function of the nucleons is shown in Fig. [S] The scattering amplitude is proportional to 

5+A'(p)5 ^ 



Extension to the vector dibaryons is straightforward. In general case, the iS-matrix can be written as follows: 



C. P matrix. Dibaryons as P-matrix poles 

The experimental searches of exotic MSQ did not give conclusive results. Jaffe and Low [s^l, therefore, proposed 
to identify exotic MSQ with so-called "primitives", i.e., the poles of the P matrix, not S matrix. The primitives 
manifest themselves in the form of zeros of D functions on the unitary cut. Analysis of the phase shifts of meson- 
meson scattering and nucleon-nucleon scattering [s^ showed the presence of primitives with masses close to those 
expected in the MIT bag model. 

The method of Jaffe and Low is basically heuristic and can hardly be extended to modified external conditions 
e.g. at finite temperature or density. In this sense, the dynamical P matrix developed by Simonov [49*1 is perfectly 
suited for studying collective properties of nuclear matter. Earlier, it has been applied successfully for description of 
nucleon-nucleon interaction in the vacuum [i^, IsH - ls^ l . 

Since the P-matrix formalism is phenomenologically successful, it is desirable for any model to have a link with it. 
As noticed by Simonov ^49i] . the P matrix has poles at primitive masses provided the interaction between nucleons 
and dibaryons is restricted to the surface of a radius h. In the coordinate space, the form factor entering the NND 
vertex represents a delta- function 5{r — b). 

In the momentum representation, we have to use therefore, in the partial L-wave, 

^"-^'^ ^ [^) — m-^ — ■ (^"-^^ 

The partial wave scattering amplitudes of the pointlike nucleons and dibaryons discussed so far already have the 
correct threshold behavior. The denominator (kb)^ in Eq. (IVII.7P is thus introduced not to violate that feature. The 
additional dependence on s/(4m^) is introduced to simplify calculations of the dibaryon self-energy operators. The 
common factor is chosen to fulfill 

J-i(s = 4m^) ^ 1. 

In the nonrelativistic limit, J- Lis) = 1 for pointlike interaction with 6 = 0. 



(1. 




FIG. 8: Four-point Green function of nucleons. 
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Let us consider function 



where jL{z), Ubiz), and h^^\z) = jbiz) +iyL{z) are, respectively, the 1-st, 2-nd and 3-rd spherical Bessel functions 
defined as Ref. 89]. One has, for z > 0, 



jL{z)-jL{-z) = 

Vhiz) - yd-z) = 



0, even L, 

2jL{z), odd L, 

2yL{z), evenL, 
0, odd L 



and, as a consequence, 

fdz)- !l{-z)^2zjI{z). (VIL8) 
The form factor modifies the discontinuity of 11/ jp(s): 

/ S f {2L+iy.l Y fL{kb)~fL{-kb) 

- ^uH<ic,Pc)[^) [-jkEy-) • (vn.9) 

Hjjp (s) can be reconstructed using the dispersion integral. 

In order to simplify calculations, we replace the ratio fijp {qciPc)/{kb)'^^ by its threshold value. Such a procedure 
can be viewed as a redefinition of the form factor J-l{s) which is a phenomenological element of the model. One can 
easily verify that the discontinuity has the correct form for the expression 

where fL{s) = /i(p*(s)). Iljjp{s) defined this way is analytical function in the complex s-plane with the cut 
(4m^, -f oo). On the first Riemann sheet, 11/ jp (s) ^ ior s oo. Equation for 11/ jp (s) can be written as 

«//pn//p(s) = (vii.io) 

with coefficients kjjp listed in Table IVlIl 

P matrix is defined as the logarithmic derivative of two-nucleon radial wave function: 



p_ d\nxL{r) 



dim 



(VII.ll) 

r—b 



Here, we suppress indices of the channel except for L. 

One can express S matrix in terms of P matrix and A(s), P matrix in terms of S matrix and A(s), and A(s) in 
terms of S and P matrices. 

Firstly, using Eq. (jVILlip and Eqs. (|VII.6P and (|VIL10p one gets 

S ^ -"D^'^-'^'D'"' (VII.12) 



kA — izjL{z)h^^\z) 
Kk + izjL{z)h^l\z) 



(VII. 13) 



where z = kb (see e.g. |9(|, problem 82). 
Secondly, P matrix can be found to be 



^ ^ {zh'^;^\z)Y + s{zh^l\z)y 

h^^\z) + Sh^j^\z) 

= Pfree+K-^I^~\ (VII.15) 
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where 

Pfree = — ' (VII.16) 

Equation (jVII.lSP is equivalent to Eq. (jIV.27P with = 0. P matrix can be modified in general by adding to P 
matrix of the noninteracting particles new poles with positive residues. The negative ones are not allowed since they 
bring the corresponding D function out of the class of generalized R functions. The free P matrix has residues 



aL 



2 



J' 

where gq = it, ai = 4.49 etc. are zeros of the Bessel functions. 

Equation (jVILlSp can be used for the extraction from experimental phase shifts the information on masses of the 
compound states which show up as P-matrix poles on the real axis of the complex s-plane. 

Thirdly, A is expressed through S and P matrices: 

kA = -izjl{z) s_i (VII.17) 

= {P-PfreeV^- (VII.18) 

Poles of A are the CDD poles. 



D. P matrix and Castillejo-Dalitz-Dyson poles 

Castillejo, Dalitz and Dyson [85] introduced the so-called CDD poles to parameterize the uncertainty of the am- 
plitudes which satisfy the correct analytic properties in the complex energy plane, the unitarity relation and are 
solutions of the scattering Low's equation. The physical meaning of the CDD poles has been clarified by Dyson [s^ 
who showed, using one of the earlier versions of the Lee model [83, that CDD poles correspond to bound states and 
resonances. New type of CDD pole is associated with primitives [93 ]. 

The CDD poles are poles of A and D functions defied by the denominator of Eq. (|VIL6p : 

D{s) ^ K{s) -n{s), (VIL19) 

where the indices {I,J^) are suppressed. The CDD poles are localized between zeros of A(s), i.e., between M^jjp 
and M^^^jjp. 

D function we discuss is the generalized R function 85]. It has no zeros on the first Riemann sheet of the complex 
s-plane. It has no zeros and on the real half-axis (— cx),so), corresponding to bound states, also provided -D(sq) < 
and So < A/^/jp- 

Resonances are simple roots of equation 

D{s) = 0, (VII.20) 

located on the second Riemann sheet below the unitary cut. In the limit of small ga, roots of Eq. (|VII.20p are localized 
in the neighborhood of s = M^jjp . 3?s gives the resonance mass, 3s determines the decay width = g'^'^D{M^)/Ma- 
At the CDD poles S{s) = mod (tt), the phase shift crossed the level with a positive slope. Indeed, let be a 
CDD pole. Equations (|VII.5|) and (jVII.6|) imply A~^(s-y) = and A~^(s-y)' < 0. Expanding the D function around 
s — and using Eq. (jVII.6[) . one finds 

6{s-yy = -^D{s^)A-\s^y > 0. 

The same behavior is implied by the Breit-Wigner formula, according to which on isolated resonances phase shift 
increases by tt. In the theory of potential scattering, growth of phase is associated with attraction. 

Models such as the Dyson-Lee model [s^ [s^l apply directly to systems where the phase shifts increase with energy. 

Nuclcon-nucleon interaction, on the contrary, points toward the dominance of repulsion, that results in decrease 
of phase shifts with increasing energy. It may seem that the Dyson-Lee models are not applicable to the problem 
of nucleon-nucleon scattering. However, this is not true. A delicate generalization allows one to extend the CDD 
methodology of constructing analytical unitary amplitudes to systems with repulsion: 
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In the earlier studies it was always assumed that '^D{s) is strictly positive [85|,|86[. Simonov [49'| showed that the 
softening of this constraint to 9i?(s) > is sufficient to ensure that systems with repulsive interactions got into the 
class of Dyson-Lee models. 

The scattering of two nucleons through compound 6QS, dibaryons, can be described with a form factor J-'l{s) that 
has simple zero(s) at s = Sp > sq = 4m^ like in Eq. (j VII. 71) . As a consequence, ^D{s) ~ (s — Sp)^ in the vicinity of 

Sp . 

If ^D{sp) 0, the phase touches at s = Sp one of the levels 5{s) = mod (tt) without crossing. 

However, if the real and imaginary parts of the equation (|VII.20p vanish simultaneously at one point, the phase 
crosses one of the levels with negative slope. One can verify this by expanding the D function in the neighborhood 
of s = Sp. Given the equation (|VII.6p and the conditions 5£'(sp)" > and diD{sp)' > 0, we obtain the following 
inequality 

In the theory of potential scattering, negative slope of phase is associated with repulsion. 

The compound MQS show up as poles of the P matrix. By their physical characteristics, poles of P matrix can be 
divided into two groups: 

The first group is a group of bound states and resonances: 

Provided that D{sq) > there is always one bound state. Other bound states are generated by the compound 
states with masses Ma < ^/sq. 

A distinctive feature of resonance is the condition J-l{s) 7^ in the neighborhood of s = AI^j jp on the second sheet 
of the Riemann surface. Equation (|VII.20p can then be used to find the simple poles of the scattering matrix. 

The roots of Eq. (jVII.20p . which lie on the real half-axis (— oo,so) of the second sheet of the Riemann surface, 
define virtual states. Their position is determined by the masses of the compound states and interaction. 

The poles of the second group are associated with the roots of Eq. (|VII.20p on the unitary cut in the neighborhood 
of s = M^jjp . Such poles do not lead to resonant behavior of the scattering cross sections, they do not generate poles 
in the S'-matrix and are not resonances. Jaffe and Low have proposed to call such states "primitives" [l^. When 
resonance hits the unitary cut, its singular effect on the S matrix disappears. In contrast to the resonances, primitives 
drive phase shifts down and mimic repulsion. 

Resonances and primitives do not exist as asymptotic states. In the diagram technique of Feynman, propagators of 
primitives are multiplied by form factors of the vertices. Such a combination has no poles at s = M^jjp . Primitives, 
thus, do not propagate, although produce a direct impact on the dynamics. 

In the general case, therefore, CDD poles correspond to bound states, resonances, and primitives. The last ones 
are P-matrix poles associated with the zeros of D functions on the unitary cut, which do not manifest themselves as 
poles of S matrix. The primitive-type CDD poles arise in systems with repulsion. 



E. Nucleon-nucleon scattering in the vacuum 



Techniques developed in Ref. 85] can be used to parameterize the nucleon-nucleon scattering amplitudes by func- 
tions with correct analytical properties within the QCB model. 



1. ^Si channel 

The most striking feature of the '^Si nucleon-nucleon phase shift is the vanishing of the phase shift at Tiab — 354 
MeV. S'-matrix according to Eq. (jVII.6p is unit in two cases: A(s) = 00 (a CDD pole) or 3n(s) — 0. 

At the CDD poles S{s) — mod (tt) with a positive slope, whereas the slope at Tiab — 354 MeV is negative. The 
CDD poles generate resonances for 5n(s) < 0. 

Let us consider the primitive-type CDD poles. The imaginary part of n(s) is related to dibaryon decay width, so 
it must be nonpositive. According to Eq. ljVII.lOp . 

QkU{s) = -zjLizf < 0. (VII.21) 

5Kn(s) has second order zeros, connected to the simple zeros of jhiz). If A(s) — 3?n(s) is finite, then S{s) = 
mod (tt) but the phase does not cross the horizontal axes 0, ±7r, .... If, however, A(s) — 3fjn(s) = and s is a simple 
zero, the phase crosses the horizontal axes 6{s) ~ mod (tt) with a negative slope. 
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A negative slope of the phase is associated with the repulsion. If we want to describe the repulsion, we have to 
fulfill Eq. (|VII.20p both for real and imaginary parts at real s — Sp > Sq = 4m^. The condition 3n(s) < gives third 
equation 

5n(sp)' = 0, 

which is, however, fulfilled automatically. 

Tiab = 354 MeV corresponds to fc = 408 MeV. Equation '^D{s) = simplifies to kb = n. We thus fix 6 = 1.52 fm. 
This value is close to 6 = 1.41 — 1.44 fm obtained in Ref. (49|]. The zero of the phase alone thus determines n(s). 

In the S'-wave, 

sin 2; ■ 
KU{s + iO) = — e*^ 

and we observe that 5Rn(s) — for s ~ Sp. Eq. (jVII.20[) simplifies therefore to A(s) — 0. The inverse A function has 
a pole at s = Sp. Equation A(s) — allows to find the dibaryon mass M — 2\/W + rn^ — 2047 MeV. This value is 
close to the dibaryon mass M — 2071 — 2081 MeV found in Ref. [i^. It is related to the primitive-type CDD pole at 
M = 3203 MeV as shown below. 
We set 

7 = -nnUis = Sd) > 0, 
where Sd = is the deuteron pole. In the case of one dibaryon, one can write 

^-lA-i = cp{^ ^) - (VII.22) 

s- Sp Sd~ Sp 7 

where Cp is a free positive definite parameter such that KCpVp — g . 

The constant term entering A~^ describes the four-fermion contact interaction between the nucleons. 

The parameterization (|VII.22p ensures the existence of zero of D function and of a pole of S matrix at s = s^, 
respectively. 

To eliminate unphysical zeros of the D function, one has to constrain Cp. We write the D function in form of the 
dispersion integral 

D{s) = A(s) - - / ^-^ds'. 

TT 74^2 S' - S-l\) 

The imaginary part of D can be represented as follows 

^ 1 /•+°°5II(sO_, , 

Since Qn(s') < 0, the bracket is positive for Cp > 0, in which case D{s) has no complex roots in the complex s-plane 
modulo the half-axis (— cxi,so)- 

It remains to check the half-axis {—00, so)- The derivative 

A'js) 1 -sUjs') ^ , 
^ 7 Tl^Cprp / — —ds 

(S - Sp)2 TT J,^ (s' - S)2 

is positive below the threshold. D{s) crosses the real axis sA s = Sd < sq. This is the unique root of the D function 
provided A(s) has no poles at s < sq- Let us investigate therefore roots of A~^(s). Since K~^{sd) = — 1/7 < 0, 
A^^(s)' < and A^^(s) has no poles dX s < Sp by the construction, the condition A^^(— 00) < is sufficient to 
exclude unphysical roots below the threshold. The parameter Cp should finally be subjected to the constraint 

< C„ < C-'^ = ^_V^^ (yjj 23) 

^ ITp 

In Fig. [9] (a) we show the '^Si phase shift versus kinetic energy of the proton in the laboratory system for Cp = 
0.9c™^'^. In Fig. ini (c) ^D{s) and 5-D(s) as functions of T are shown. We recall that for pn system s = sq + 2mnTiab 
where m„ is the neutron masses. 3?I?(s) has one zero below the two-nucleon threshold, which corresponds to the 
deuteron. 

If the number of dibaryons is greater than one, the constraint IVII.231 reads 

p Sp-Sd 7 

where Cp are positive. 
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FIG. 9: (color online) '^Si and ^So scattering phase shifts in radians (upper panel) and real and imaginary parts of the D 
functions (lower panel) versus the proton kinetic energy. The long-dashed curves are parameterizations of the QCB model from 
Ref. [4^, the solid curves are parameterizations of the relativistic QCB model. The experimental phase shifts [93] are shown 
by the circles. 



2. ^ So channel 



The scattering phase vanishes at Tiab = 265 MeV. Arguments similar to the ones used for the ^Si channel lead us 
to the estimates h — 1.76 fm and M = 2006 MeV. The numbers are close to those obtained in Ref. [i^ h = 1.44 fm 
and M = 2112 MeV. The CDD pole is localiz ed at M = 2916 MeV. 

The inverse A function has the form of Eq. (jVII.22l) with Sd replaced by sq = 4m^, where m is the average mass of 
proton and neutron. Near two-nucleon threshold, kD{s) — —7 + 1 + ikb + ... From other hand, D(s) ~ 1 — i6{k) + ... = 
1 — ika + where a = 23.56 fm is the scattering length. One has to require 

7 = -Kn(so) + ^, (VII.24) 

with n(so) = -1. 

D function of the ^So channel has no roots for the complex values of s and its derivative is positive on the real 
half-axis below two-nucleon threshold. In order to avoid the unphysical poles of S matrix, it is sufficient to require 

^^(-oo) = cpr^ 1 < nD{so) = 1 - 7 < 0. 
Since 7 > 1, the second inequality is fulfilled. The first one gives 

Cp < min(c--^ (VII.25) 
7 — i 

where dp'^^ — {sp — so)/{'jrp). Since b <^ a, it reduces to Cp < c™^'^. 
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In Fig. [9] (b) we show out fit of tlie ^Sq phase shift versus the proton kinetic energy with Cp = 0.9c™'^'^. In Fig. 
ini (d), real and imaginary parts of the D function are shown on the real Tj^b-axis. 5ftZ3(s) has no zeros below the 
threshold. ^D{s) is positive definite above two-nucleon threshold. At s = Sp both the real and imaginary parts of 
the D function vanish. 

3. "^Pi and ^ Pi channels 

The ^Pi and ^Pi phase shifts dot not cross the real axis. So, we have to use other methods to estimate the parameter 
b and the dibaryon mass. 
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FIG. 10: (color online) ^Pi and ^Po scattering phase shifts in radians (upper panel) and real and imaginary parts of the D 
functions (lower panel) versus the proton kinetic energy. The solid curves are parameterizations within the relativistic QCB 
model. The experimental phase shifts [o^l are shown by circles. 

The inverse A function can be taken in the form of Eq. (jVII.22p with replaced by sq = Am? . At the threshold, 
kD{s) — —7 — Kn(so) + ^k^b^ + .... From other hand, D{s) ^ 1 — iS{k) + ... = 1 — ik^a + where a = —1.5 x 10^'' 
MeV~^ for the ^Pi channel and a = —3 x 10""^ MeV""^ for the ^Pi channel. To provide correct threshold behavior 
one has to require 

7 = -^U{so) + ^. (VII.26) 
Since a is negative, kD^sq) = —kII(.sq) — 7 > 0. In order to remove unphysical poles, it is sufficient to require 

< kD{-00) = cprp'^ 1 < kD{so) = -Kn(so) - 7. 

Sp — Sq 7 

In the case of 7 < 0, this equation provides no constraints on Cp. In the case < 7 < —kII{so), k~^A~^{sq) = — i < 0, 
K^^A^^(— 00) = -j^^z^ — ^ > 0, and A~^(s)' < 0. There is therefore a pole in the D function, and zero of the D 
function on the real half-axis (—00, sq). The case of < 7 < —kII{so) is thus not useful. 
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Let us assume 7 = 0. The value of b for the ^Pi channel can be found to be b = 1.51 fm. The dibaryon of mass 
M = 2214 MeV that makes A(sp) to vanish ensures the crossing of the level S{s) = — tt at Tiab — 732 MeV. In such a 
case, the repulsion is largely overestimated. The linear interpolation of the experimental data e.g. gives S{sp) ~ — 1. 
If M 2214 MeV, the phase remains negative and touches one of the levels S{s) — mod (tt) at s = Sp. A better 
description is possible for lower value of b and greater value of M, respectively. However, lower values of b require 
positive 7 which is not acceptable. 
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FIG. 11: (color online) '^Pq scattering phase shifts in radians (upper panel) and real and imaginary parts of the D functions 
(lower panel) versus the proton kinetic energy. The solid curves are parameterizations within the relativistic QCB model. The 
experimental phase shifts fo^ are shown by circles. 



We thus consider more general expression for the imaginary part of n(s): 



(VII.27) 



where x is a, parameter. n(s) is restored from the dispersion integral. The explicit analytical expression for n(s) can 
be found e.g. with the help of MAPLE. For x = 00, we recover previous formulas. A finite x leads to reduced values 
of n(so) in the absolute value and the parameter 6, respectively, keeping 7 negative. For L — 1, 

KHiso) = -^(-3 + ix"^ -2x^ + 3(1 + xf e"^^). 

One can notice that < — Kn(so) < 1/3. 

The same arguments as before lead us to the conclusion of absence of unphysical poles of the S'-matrix in the 
complex s-plane for negative 7. 

A reasonable fit in the "^Pi channel can be obtained without the CDD poles and primitives. Cp = implies that 
the only important interaction is the contact four-fermion interaction described by 7. We use 7 — —1/20, b — 0.93 
fm and x'^ = 1/3 which corresponds to kII{sq) — —0.028. 

In Fig. [TO] (a) we show phase shift of the ^Pi channel. Shown in Fig. [TO] (c) are real and imaginary parts of the D 
function. 
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The phase shift in the ^Pi channel behaves similar to the phase shift in the ^Pi channel. We use therefore precisely 
the same set of the parameters. The scattering length appears to be a factor of two lower that the experimental one, 
however, the overall description of the data is better than with the increased scattering length. In Fig. [10] (b) the 
phase shift in the ^Pi channel is shown. In Fig. [Tn](d) the real and imaginary parts of the D function are shown. 



4- ^Pq channel 

The parameterization of the phase shift in the ■^Pq channel is similar to that in the ^Sq channel. Here, we do not 
modify 5n(s) and set xr = oo. The values of 6 = 3 fm and the dibaryon mass M = 1969 MeV are found from the 
position Tiab = 186 MeV of zero of the phase shift. The primitive-type CDD pole appears at M = 2650 MeV. The 
parameter (cf. Eq. (jVII.26P ) 

can be fixed by the scattering length a — 2.8 x 10^'' MeV"'^. We consider, however, the scattering length as a fitting 
parameter. The general agreement with the data is better when slope of the phase shift at Tiab = is two times lower 
than the experimental one. Fig. [TT] (a) shows the phase shift as compared to the experimental data as functions of 
the proton kinetic energy. Shown in Fig. [TT] (b) are the real and imaginary parts of the D function. 



F. Discussion 



In Sects. V-VII we constructed and analyzed the effective Lagrangian to describe the two-nucleon forces based on 
their s-channel dibaryon exchange mechanism. Dibaryons in our approach are the primitives, and not resonances. 
They manifest themselves as poles of P matrix or, equivalently, as zeros of D functions on the unitary cut. Resonances 
with the energy-dependent width have similar properties, when the width vanishes at the resonance mass, while far 
from the mass of the resonance the width is finite. The effective Lagrangian was formulated for the nucleons, which 
are well-defined asymptotic states, and for the primitives, which are decoupled of the two-nucleon channel on the 
mass shell, and, apparently, are not asymptotic states. Note that although the resonances are not asymptotic states, 
they are often included in effective Lagrangians. Despite the unusual properties of the primitives, there is no intrinsic 
difficulty in applying the diagram technique for calculation of various processes involving primitives, based on their 
effective Lagrangian. 

In Sect. V we gave a complete list of relativistic two-nucleon currents with the angular momentum J = 0, 1, through 
which the primitives are coupled to two nucleons in the low- J partial waves. We found five primitives with the lowest 
quantum numbers (/, J^) = (0, 1^), (1, 0^), and (1, 1~) and discussed their effect on the nucleon-nucleon interaction 
in Sect. VII. In the various versions of the QCB model, the 5'-wave nucleon-nucleon phase shifts are fitted by the 
primitives excellently and we confirm it again. We considered the P-wave scattering also. We found good qualitative 
agreement with the experiment. The description of the channels '^Pi and ^Pi requires the contact four-fermion 
interaction only, without primitives, whereas the channel ^Pq gives an evidence for the primitive. We suppose more 
accurate parameterization of the experimental phases can be achieved by means of the s-channel dibaryon exchange 
only without invoking the t-channel meson-exchange mechanism. 

P-matrix analysis has been successful in describing the various processes, including nucleon-nucleon interaction. 
This formalism is built into the effective Lagrangian with an appropriate choice of nucleon-nucleon-dibaryon form 
factors and four-fermion vertices. 

In the limit of low density we reproduce results of optical potential model. Lorentz structure of the nucleon self- 
energy operator is different from the OBE RMF models. In our approach, vector and scalar parts of the self-energy 
operator are of the same order of magnitude. 

Solutions of the self-consistent system of equations for the nucleon and dibaryon propagators and G-matrix in 
nuclear matter will be presented elsewhere [9^ . 



VIII. MODELS OF NUCLEAR MATTER 



The nuclei are known to be stable under normal conditions on Earth, while the collection of nuclei forms the ground 
state in QCD at finite baryon charge. 
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The range of densities and temperatures p ^ 10'^^ cm~^ and T ^ 10^^ K, where the phase transition in nuclear 
matter takes place with chiral symmetry restoration and deconfinement, can be studied experimentally in heavy- ion 
collisions. Supercold nuclear matter exists in the interiors of massive neutron stars. 

Quark matter appears in compression of nuclear matter. New forms of nuclear matter can exist at intermediate 
densities, such as hyperon matter, crystalline neutron matter, pion and kaon condensation and a Bose condensation 
of dibaryons. 

In the field of atomic physics, modifications of the properties of ordinary atoms at density p ~ 10^" — 10^^ cm~^ 
and temperature T ^ 10** K are studied since the mid-1930's (for a review see [96l, |97|). 

To ensure the stability of nuclear matter under conditions T < 10 ^ 20 MeV and p ~ 0.16 0.30 fm~'^, the NN 
forces must satisfy necessary and sufficient saturation conditions (see, e.g., the review of F. Calodgero and Yu. A. 
Simonov jiSj ) . For local and velocity-independent potentials these conditions require a strong repulsive core at small 
distances, r < 0.4 fm. 

The quark-hadron phase transition in nuclear matter has been discussed first by Zel'dovich (99| . By comparing 
pressures of the ideal Fermi gases of nucleons and quarks, he arrived at a conclusion that the composite structure 
of the nucleons requires a short-range repulsion between the nucleons. His discussion was based on the Le-Shatelie 
- Brown principle. Earlier he showed that the stiffest, consistent with causality nuclear matter EoS corresponds to 
exchange by a massive vector particle between the nucleons (lOOj |. The w-meson exchange mechanism is the most 
typical and challenging feature of the modern OBE models. 

As discussed above, at short distances NN force becomes nonlocal - in OBE picture due to finite size of nucleons 
and mesons - and in the QCB picture due to formation of MQS. It is remarkable that nonlocal NN interaction, even 
when everywhere attractive, easily passes the requirements of saturation, while in the local OBE case one needs a 
very large coupling of w-meson, producing repulsive core. 

The realistic OBE models of nuclear matter split into three groups: 

Relativistic mean field (RMF) models are based on the mean-field approximation to account for the meson-exchange 
interactions between nucleons, a- and w-mesons ensure the saturation property. The mean field models are the 
phenomenological ones. 

Dirac-Brueckner-Hartree-Fock (DBHF) approach is a field-theoretic scheme, which is based on accurate experimen- 
tal data of nucleon-nucleon interaction. Nuclear many-body problem is solved microscopically by the summation of 
wide subclass of Feynman diagrams. The scheme claims to be parameter-free description of the nuclear matter EoS. 

The variational approach is proposed in (26j . This uses the potentials of the nucleon-nucleon interaction, whose 
characteristics are extracted from the nucleon-nucleon scattering and descriptions of light nuclei. Ground-state wave 
function of infinite nuclear matter is constructed with the use of the variational principle. Nuclear EoS is claimed to 
be parameter-free as well. 

The DBHF models and the variational approach are successful in the description of nuclear matter at saturation 
starting essentially from the first principles, however, with an additional tuning, e.g., of 3iV forces in the variational 
approach. 

In addition to the above three methods, the non-relativistic BHF calculations [lOl], phenomenological theories with 
density dependent interactions such as Gogny or Skyrme forces ^102,] and various modifications of the RMF scheme 
are widely used for modeling nuclear matter also. 

A critical test of nuclear models is possible on the basis of astrophysical observations on the properties of massive 
neutron stars. 

The most massive pulsars reported in the lit erat ure are PSR B1516-I-02B with the mass of 1.96^g'52 Mq and PSR 
J1748-2021B with the mass of 2.74 ± 0.22 Mq jlOSj . The mass of rapidly rotating neutron star in the low mass X-ray 
binary 4U 1636-536 is found to be M = 2.0 ± 0.1 Mq Q. The mass of the X-ray source EXO 0748-676 is limited 
to M > 2.10 ± 0.28 Mq. These observations indicate clearly that the /3-equilibrated nuclear matter EoS is stiff and, 
furthermore, eliminate the soft EoS constructed, e.g., on the basis of the classical Reid soft core model |26|]. 

The new degrees of freedom lower the energy of matter and make EoS softer. Scenarios with phase transitions shift 
the masses to the region of even smaller values. Existence of quark matter and other forms of exotic matter in the 
cores of neutron stars was questioned [l^ (see, however, (33|). 

The so-called realistic models of neutro n sta rs usually ignore hyperon channels, for example, the reactions — > 
n + e + De- In the RMF models [13, Il04| . which take hypernuclear data into account, the inclusion of the (3- 
equilibrium balance for hyperons reduces the limiting mass of neutron stars by 0.5 0.8 Mq. These results are in 
agreement with recent calculations [29l - [3l| . The /^-equilibrium brings thus difficulties to reproduce the observed high 
masses of neutron stars. 

The problem can be solved by assuming the existence of light weakly interacting bosons, the possibility of the 
existence of which is discussed in some generalizations of the Standard model [32^ . A more conventional and perhaps 
more realistic approach is to modify the baryon-baryon interaction at short distances, taking the quark-gluon degrees 
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of freedom into account. It is based on the formalism of the P matrix and QCB model discussed in the previous 
sections. 



A. Relativistic mean field models 

Development of many-body theory in nuclear physics was motivated by the problem of describing heavy nuclei 
and, since 1970's, neutron stars. The change in mass of the nucleon in nuclear matter has been introduced in the 
pioneering work of Walecka |24|. Mean- field (MF) and relativistic mean field (RMF) approximations were proposed to 
describe the dense nuclear matter self-consistently to all orders in density, but in tree level of perturbation theory. MF 
and RMF were very successful in describing the basic properties of finite nuclei and nuclear matter near saturation 
density. RMF models are currently the most popular models of nuclear matter, although they do not provide a strict 
connection with the phenomenology of nucleon-nucleon scattering. 

Since mid-1970's, RMF models have been improved in many directions: by including nonlinear self-interactions of 
a-, p-, and (5-mesons, introduction of the density-dependent coupling constants and vertices with higher derivatives, in 
order to describe the energy, radius, surface thickness of nuclei, neutron skins, spin-orbit interaction of nucleons, as well 
as restrictions on flows obtained in heavy-ion collisions, and fulfill constraints imposed by astrophysical observations 
of neutron stars. 

In RMF models a significant decrease of the effective mass of the nucleon was found at a density of saturation. 
This effect was predicted purely phenomenologically. The development of QCD confirms the prediction and the effect 
has attracted much attention, since the decrease of the nucleon mass is connected to the partial restoration of chiral 
symmetry. 

In the simplest version, the effective Lagrangian contains the nucleon fields and a- and mesons: 

- ^^tJ.y^f^iy + ^mlujf,UJ^„ (VIII.l) 

where J"^,. = V,.w^ - V^w^, V = 7^V^, etc. 

Property of saturation is easily seen from the following argument: Suppose that low-density attraction of the scalar 
field dominates the uj exchange. With increasing density, the scalar charge of the nucleons, generating a scalar field, 
is suppressed by the Lorentz factor (The scalar charge of a particle placed in a box disappears in the ultrarelativistic 
limit. The vector charge, of course, remains constant.) In infinite-density limit, cr-meson decouples from the system of 
nucleons, but vector fields survive. As a result, the attraction at low densities is replaced by the w-induced repulsion 
at high densities. In some intermediate-density regime, repulsion balances the attraction, this density is the density of 
saturation where the substance is in equilibrium. The coupling constants go-, guj, and the cr-meson mass are determined 
from fitting the properties of nuclear matter at saturation density. 

RMF models offer low binding energy of nuclear matter at saturation due to cancellation between the strong 
repulsive vector and strong attractive scalar potentials. These two potentials are equal in magnitude to about a 
quarter of the nucleon rest mass. The scalar part reduces the nucleon mass and increases the infiuence of small 
components of the Dirac bispinors. 

RMF models allow the most precise description of the nuclear properties with 5-15 parameters. 

B. Dirac-Brueckner-Hartree-Fock approximation 

DBHF approximation scheme is based on realistic meson-exchange models. As compared with the RMF model, it 
takes additional requirements of fermionic statistics and in-medium modification of two-body T-matrix into account. 
DBHF successfully reproduces the properties of nuclear matter at saturation mainly from first principles, using the 
interaction, based on data from nucleon-nucleon scattering. 

Shown in Fig. [T^] is the Bethe-Salpeter equation for the in-medium two-body scattering T-matrix (called reaction 
G-matrix). The difference from the vacuum equations is associated with the change of nucleon propagators in the 
loop and Pauli blocking, which restricts the admissible regions in the phase space of nucleons. 

The nucleon self-energy can be found by solving equation shown in Fig. 1131 

Equations in Figs. [TH and [T3] constitute selfconsistent system of equations of the DBHF approximation. 
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FIG. 12: System of equations for in-medium T-matrix. V is meson-exchange potential. The bold lines in the loop are in-medium 
nucleon propagators. 




FIG. 13: System of equations for in-medium nucleon self-energy operator. The bold lines in the loops are in-medium nucleon 
propagators. T is in-medium T-matrix of Fig. 1121 

Examples of the DBHF calculations can be found in Refs. DBHF scheme is phenomenological only 

marginally: The ladder approximation is used to solve the relativistic Bethe-Salpeter equation, negative-energy nu- 
cleons are neglected, nucleon resonances are neglected also. 



C. Variational approach 

Estimates of G-matrix with realistic potentials show that G-matrix has a sharp peak with respect to the difference 
of distances between the nucleons in the initial and final states. This means that in the process of scattering the 
distance between nucleons is not much changed, and for the variational calculations the dependence of the trial wave 
function on the relative momentum of the nucleon pairs can be neglected. The trial wave function has the form 

i<j i k 

where is projection operator to state of the pair ij with the orbital momentum A is the antisymmetrization 
operator, (j>{rk) are uncorrelated wave functions of the nucleons. The correlation function /^(r) is equal to one when 
the distance between the nucleons increases, but less than unity, when r is reduced. 

The Jastrow wave function (jVIII.2l) is used for the calculation of the ground state energy from the cluster expansion 
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etc. Here, H2 is two-body Hamiltonian, is two-body wave function. 

The smallness of the correlation effect in the two-body nucleon wave function, 



Po 



^ J dr{l - f{r)f ^ 0.1, (VIII.5) 



at the saturation density po = 0.16 fm ^ makes it possible to truncate t he cluster expansion. 
Recent examples of the variational calculations can be found in Refs. jl09l . Ilici| . 
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D. In- medium modifications of hadrons 



In-medium modification of tfie properties of hadrons affects tfie nuclear matter EoS. In the RMF models, vector 
mesons play a special role. The interaction mediated by vector mesons at high densities leads asymptotically to the 
stiffest EoS of nuclear matter. The effect of the vector meson on the energy density of nuclear matter can be estimated 
by averaging the Yukawa potential jlOCll | 

£/ = dxidX2p(xi)|^^p(x2), (VIII.6) 

where p(xi) = p(xi) = p is baryon number density, r = |x2 — xi|, g is coupling constant, fi is vector meson mass, and 
V is normalization volume. A simple integration gives 

= (VIII.7) 

The RMF models are used to study the modifications of the properties of vector mesons in nuclear medium 
Ill3j |. Assuming that the mass of mesons depends on the density, the contribution of the w-meson to pressure can be 
found to be 

p,=^ri-^I^V (VIII.8) 



2^2 \^ fJ- dp 

This equation shows that the pressure decreases if the meson mass increases with density. The problem of modifying the 
mass of vector mesons in nuclear matter was studied in experiments on the dilepton production in heavy-ion collisions. 
Simulations of the dilepton production using transport models and the subsequent comparison with experiment gives 
a clear indicati on o f the presence of collisional width, however, a noticeable shift of the masses of p- and w-mesons is 
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IX. NUCLEAR MATTER WITH MULTIQUARK BARYONS. NONRELATIVISTIC TREATMENT 

For two-nucleon contrib utio n to the nuclear matter one starts with the Bethe-Brueckner series which is summed up 
by the operator equation [llSj 

G2{W) = V2-V2Q2e2^G2{W) (IX.l) 

where G2{W) is the reaction matrix of two nucleons in nuclear matter V2 is the two-body interaction operator 
(potential), Q2 is the projection operator on states outside of the Fermi sphere, and finally 62 is energy denominator 

e2|pn >= {Ep + En - W^)|pn > . (IX.2) 

Eq. (|IX.1I) is an equivalent of two-body i-matrix (Lippmann-Schwinger) equation for the case when two nucleons are 
inside nuclear matter. 

One can also write equations for the pair wave functions. For the plane wave $pn and actual w.f. ^pn one has 

*pn $pn - Q2e2^^G2$pn (IX.3) 
iKR, 

Separating out the cm. motion, ^'pn = ^-^7 — V'prij a-nd keeping only relative coordinates in V'pnj^'pn = exp(ikr), 
one has the Bethe -Goldstone equation (for nonlocal in general interaction) 

V'pn(r) = e'''^ - j F2{r, v')V2{v', r"; E)^pn{r")dr'dr" (IX.4) 



where 



^Q2(k,K)[e2(k,K)]-iexp[zk(r-r')], (IX.5) 

g2(k,K)|p,n> = I IP'">' l5K±k|>fc^, (j^g^ 
I 0, otherwise, 

ealpn) = [ii;(^ + k) + i?(| - k) - T4^]|p, n) = e2(k, K)|p, n), (IX.7) 
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and 



K = p + n, k=i(p-n). (IX.I 



From tppn one easily finds the matrix G 



l,m|G2|k,n) = l<5j^j^,(k'|G2|k) (IX.9) 



(k'|G2|k)= J e-'^''^V2{r,r')^pn{r')dvdr'. (IX.IO) 

Inserting the separable interaction V2(r,r') ~ Vhqh ~ ^"e^'-e'^ ' obtains from pX.4[) the following system 
of algebraic equations 

{Ui>) = Mk) - ^ A,,.(/,,v) (ix.ii) 

where we have used notations 

(M) = J drMr)^pr{r), Mk) ^ J dre^'' Ur) , (IX.12) 

_ J Mr)F2{r,r')U,{r')drdr' _ A^^l ..^ ^ 

In a similar way one easily calculates (k'|G|k), 

(k'|G2|k)=5:/;^M^. (IX.14) 

Keeping only one MQS level in pX.ll|) . pX.14|) one has 

(k'|G2|k) = ^j^^^g^ 

The case of more MQS levels is discussed in the next sections. Note, that E in all equations in this paper denotes 

the total energy of two-nucleon system, including its masses, which finally in (|IX.20p coincides with W, and one- 
particle energies E{p) also include rest mass in kinetic term. To make pX.ip selfconsistent, one should also define 
the single-particle energies En , Ep in (IIX.2p through the same matrix operator G2 [ llSj l . 

E{p)^T{p) + U2{p,p) (IX.17) 

f/2(p,p) = {PP'\G2{W - E{p)+E{p'))\pp')a (IX.18) 

P'<Pf 

where subscript a implies antisymmetrization of state (pp' > . 

To calculate the energy per nucleon one can use the BHF approximation 



where we have defined 



({/2)= (PP'|G2(£;(p)+i?(p'))|PP')a. (IX.20) 

p,p'<PF 
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We now turn to the contribution of higher MQS states, and start with the case, when the only interaction present 
in nuclear matter is the MQS composed of 3A^ quarks, N > 2, while all other interactions with N' ^ N are absent. 
Denoting the corresponding reaction matrix for N nucleons G n , one can write similarly to (jIX.l|) 

Gn{W) = - VNQNeJj'GNiW) (IX.21) 

where Qn, sn are natural generalizations of Q2, 62 and Vn is 

^-^ E — El, 

V 

It is clear, that generalization of Eqs. (jIX.lHIIX.16p to the case of iV-nucleon interaction is straightforward if one 
replaces there one relative momentum k by the set of — 1 relative momenta {kct},a = 1, ...TV — 1, and similarly 
generalize (32(k, K) and e2(k, K) to QAr({kQ},K), eAr({kQ,}, K). For the function Fj.^ one defines — 1 relative 
coordinates {r^jja = 1, ...N — 1, and one has instead of pX.sp 

»JV-1 N-\ 

F„({rJ,{r'J) = / J] — ^Qw({k4,K)e^i({k4,K)exp(z ^ k^r^ ~ r^)). (IX.23) 

The answer (jIX.16p for the case of one MQS made of 3A^ quarks is 

{{k'J\G^\{k,}} = /--({kLD/^dM) (JX.24) 

E — Ey + Kuu 

and A^, is expressed through F^ in the same way as in (jlX.lSp with the obtains replacement r {ra},r' {r^}- 
Finally in (|IX.19P one replaces by where (C/jv) is 



(;:^^>= E ({PJ|G^(^^)l{PJ)a (IX.25) 

{Pi}<PF 

and the set of one-particle momenta {p,}, i = 1, ---N is inside the Fermi sphere, while Wn — X^iLi ^(Pi)- 

X. G MATRIX 

Diagrams in Figs. El [71 and |8] show a closed system of equations that can be solved e.g. by iterations starting from 
the vacuum solution. 

The projection operator Q2 on states outside of the Fermi sphere influences scattering of the nucleons as described 
in Sect. IV. In the relativistic theory, the dibaryon self-energy Fig. jSlis represented in the form of dispersion integral. 
The phase space of two nucleons is restricted by Q2 and weight of the dispersion integral modifies, accordingly. The 
invariant part of the dibaryon self-energy operator can be written in the form 

iiij.{s,v) = - Bf{s',v) ' X.l 

Jso S' - S- lO 

where Bp is the Pauli blocking factor and v is the dibaryon velocity. In terms of the nucleon momenta, the imaginary 
part of 11/ jp(s) is the same as in the vacuum. The only modification of Iljjp(s) comes from Bp ^ 1 and the 
in-medium nucleon mass m* = m + 'Ss{p+) where p+ = (•\/m*^ + p^, p). 

A. Pauli blocking 

One can give explicit relativistic form of the blocking factor pX.6p . It can be used with Eq. (|X.1|) for evaluation of 
the in-medium dibaryon self-energy. 

In the rest frame, K, of the substance the Fermi sphere is described by equation 

pI+pI+pI^pI- (X.2) 
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FIG. 14: Possible configurations of the boost-deformed Fermi sphere and the decay sphere in the rest frame of two nucleons. 
Upper panel pf < p*'- Condition —p* < —y{pF + vEf) is fulfilled in the case (a). The Fermi sphere is inside of the decay sphere. 
The decays are allowed and Bp = 1. The spheres intersect each other in the case (b) when —^{pf + vEf) < — p* < ^{pf~vEf) 
and < Bp < 1. In the case (c), the relative velocity is high, 7(pf — vEf) < —p*, the spheres do not intersect each other and 
the decays are allowed Bf ~ 1. Lower panel p* < pf'- In the case of p* < 7(pf — vEf), the decay sphere is inside of the Fermi 
sphere and the decays are forbidden Bf ~ 0. The case (d) corresponds to intersection of the spheres for y{pF — vEf) < p* 
and E* < Ep/y where Bp = also. In the case (e), —p* < ^{pf — vEf) < p* and Ep/y < E* < Ep, the decays are possible 
and < Bp < 1. In the case (f), the relative velocity is high to provide y{pp — vEp) < — p*, in which case the spheres do not 
intersect and the decays are allowed Bp = 1. 



In the co-moving frame, if', of the dibaryon the decay products have momenta 

P'^+P'^+P'I=P*^- (X.3) 

We assume that the dibaryon moves in the positive direction of the x-axis. Using boost along a:, one can write equation 
for the Fermi sphere in K': 

where u > and = l{p'x + Py — P'yi Pz — P'z- Since p' belongs to both the Fermi sphere and the decay sphere 
(OOI) . E' =E* = v/m*2 + p*2 and 

y\p'^ + vE*f+p*'-p''^=pl, 

from which we derive 



, -E* ± Epy/l - 
P.± = , 



where Ep = \J m*"^ + p^p. Note that p'^_ < —p*. There exists therefore only one intersection of the spheres. The 
nucleon self-energy depends on the absolute values of the momenta, so directions of the nucleon velocities have equal 
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weights. The modified dispersion law imphes that the nucleon mass entering E* and Ep is the effective nucleon mass 

TO*. 

The critical values p'^_^ — ±p* at which the Fermi sphere touches the decay sphere give equations 

, , -e* + EfVT^ 
±p ^ . 

V 

The blocking factor is defined by fraction of the area of the decay sphere outside of the Fermi sphere. The possible 
cases are shown in Fig. [TH The blocking factor becomes 

BF = ~coae+, (X.4) 

where 

cos 9+ = max(— 1, niin(0, ^^^)). 



B. Heterophase nuclear matter with a Bose condensate of dibaryons 

Small perturbations of s-channel NN interaction, as a rule, shift primitives from the unitary cut and turn them 
into resonances. One can expect therefore that the effects of the environment transform the primitives into the 
resonances. Such dibaryons can provide deep modification of nuclear matter at densities above the critical density 
for the formation of such dibaryons on mass shell. Before proceeding to discuss realistic models it is worthwhile to 
consider a mixed dibaryon-nucleon phase in ideal gas approximation. 

Suppose that wc fill the box with neutrons, as shown in Fig. 1151 Pauli principle allows neutrons with increasing 
density consistently occupy higher energy levels. We can continue this process until the chemical potential of neutrons 
at the surface of the Fermi sphere will not be higher than half of dibaryon mass. In the approximation of an ideal 
gas, the chemical potential of fermions is equal to the Fermi energy, while the chemical potential of bosons is simply 
the mass of boson to d . 

When the Fermi energy of neutrons becomes higher than half the mass of the boson, it gets energetically advanta- 
geous for the two neutrons to merge and form a dibaryon. The critical density at which the process of the dibaryon 
condensation starts is thus determined by the mass of the lightest dibaryon. Above the critical density, the chemi- 
cal potential of the nucleons /i„ is frozen at ^i'^^^^ = tod/2 due to chemical equilibrium with respect to conversions 
NN o D. 



^ Neutrons 




FIG. 15; Illustration of transition to the heterophase state of nuclear matter with an admixture of dibaryons. 

Dibaryons are Bose particles and form a Bose condensate. Dibaryons in the condensate have zero velocity, so they 
do not collide with the boundary of the box and do not contribute to the pressure. 

On the other hand, the Fermi energy of neutrons is frozen, so the pressure does not increase with increasing density 
(see Fig. I16|) . As a result, a binary mixture of neutrons and dibaryons loses its elasticity. Nuclear matter with such 
properties can not protect neutron stars against gravitational compression and the subsequent collapse. 

Incompressibility of nuclear matter at saturation density is well known experimentally, it is certainly not equal to 
zero, so the dibaryon Bose condensate, apparently, does not exist in normal nuclei (at least in the ideal gas model and 
some simple modifications of this model, see below). Because of the charge symmetry of nuclei, /Lt„ = /ip. Assuming 
further that the shell potential of dibaryons is two times deeper than for nucleons, we can conclude that the masses 
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FIG. 16: Schematic view of EoS of the heterophase nuclear matter with an admixture of dibaryons in the absence of interactions. 



of primitives (which turn in nuclear medium to dibaryons, i.e., resonances) associated with the iVA^-channel must be 
greater than 

mo > 2eF = 1.96 GeV, (X.5) 

where sp = "/tin + AO MeV is the Fermi energy of nucleons in nuclei. The primitive masses we derived earlier from 
the NN scattering, rriD = 2047 and 2006 MeV, satisfy this constraint. 

In the nonrelativistic approximation, the evolution of the local velocity v of matter inside the newly born protoneu- 
tron star is described by Euler's equation 

at 

where $ is gravitational potential, p is density, and p is pressure. In the presence of the dibaryon Bose condensate, 
pressure remains constant, as indicated by the horizontal line in Fig. 1161 and hence Vp = 0. Gauss' law implies that 
/ dS • = 47rGM(r), where M(r) is the mass inside a sphere of radius r, so ^ 0. Thus, we conclude that 
V ^ 0. In the approximation of an ideal gas, the static solutions disappear because of formation of the dibaryon 
Bose condensate. These arguments are valid in the general relativity. Formation of the dibaryon condensate can be 
considered as a possible mechanism for the phase transition to quark matter. 

Dibaryon- nucleon and dibaryon-dibaryon interactions contribute to the pressure. EoS of the binary mixture becomes 
stiffer, ensuring the stability of neutron stars in a certain range of densities. If the formation of the dibaryon Bose 
condensate is a first-order phase transition and the density jump Ap at the phase trans ition is large enough, the 
neutron stars become unstable when the central density hits the critical value jHO, [H, Ill6j | . 

In a Bose gas of interacting bosons, a fraction of bosons is outside the condensate. These particles move, collide 
with a boundary, and contribute to the pressure. The stronger the repulsive interaction, the more stable can be a Bose 
gas of interacting bosons. With increasing density, the interaction energy is growing faster than the kinetic energy, so 
the difference between fermions and bosons becomes asymptotically less significant. Coup ling constants and masses 
of dibaryons can be constrained from the existence and stability of massive neutron stars |45| . 

In the RMF model of Ref. 45] dibaryons with masses close to the ones we derived from the phase shifts analysis 
(tod = 2047 and 2006 MeV in the and ^5*0 channels, respectively) with sufficiently weak cr-meson couplings 
(weak attraction force) appear at high densities without destabilizing neutron stars. Dibaryons with strong cr-meson 
couplings appear at densities below 0.3 fm^'^. Such dibaryons generate instability (square of sound velocity gets 
negative), providing phase transition to more dense phase of nuclear matter, e.g., quark matter. 

There is therefore the interesting connection between parameters of the lightest dibaryons and structure and stability 
of neutron stars. 



XI. NUCLEAR PHASE TRANSITION IN COLLAPSED STELLAR CORES 



The behavior of stars under a phase transition instability is determined by the jump of density A. If A exceeds the 
critical value Ac = |, the star becomes unstable at the moment the phase transition starts in the center. This result 
is remarkably independent on the specific properties of EoS. The magnitude of the density jump is the only relevant 
parameter [60j. 
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In gene ral, t he problem of stellar stability with a phase transition is more complicated as discussed, for instance, in 
Refs. [Ulllll. A star can lose stability even when A < | (subcritical regime) on conditions that the interface between 
two phases is not exactly at the stellar center and a new-phase inner core of finite size has already been formed. The 
critical size of the new-phase core depends on the stiffness of EoS which is determined by the adiabatic index averaged 
over the star 7. The lower the stability factor (7 — |) the smaller density jump is required to destabilize the star. 
The hydrostatic equilibrium of a star at the boundary of stability (7 = |) can be disturbed by phase transition with 
an arbitrary small jump of density. 



A. Equation of state 



Since our goal is to investigate only the main features of the hydrodynamic behavior of the instability we use a 
simplified approximation for EoS as described below. 
The free energy is approximated in the following way 

Fi^^ii^ E,4T,^) + AFi^,uin), n^n,d{T,0. (XI.l) 

where F is the specific (per particle) free energy, T is temperature, and n is the particle number density. The 
indices L and H correspond to low and high density phases, respectively. For simplicity, the matter is assumed to 
consist of neutrons only. The functions Fid and riid are defined by the Fermi-Dirac statistics for ideal (noninteracting) 
nonrelativistic particles and depend on the temperature and parameter ^ that failing the interaction is simply the 
chemical potential. 

The additive terms AF{n) in (|XI.ip describe the interaction and depend only on n: 



AFl,h(«) = Cl,h Ti'^^'H-i _ c^dn^/^ . (XI.2) 

Thus, Eqs. (jXI.llXI.2p define as a function of n and T parametrically (parameter ^). All the necessary for 
hydrodynamical calculations of thermodynamical quantities such as the pressure, internal energy, entropy, chemical 
potential etc. can be easily found from F by standard way. For example, the pressure P and chemical potential /i are 

/ dF\ P 



dn J n 

The region of transition between the two phases can be found from the equations of phase co-existence 

PL{T,nL) = PH{T,nH), tiL{T,nL) ^ m{T,nH), (XI.4) 

which require the continuity of pressure and chemical potentials at the phase interface. These equations determine 
the particle density of the phase transition beginning and that of its ending uh, and so the magnitude of density 
jump A = nn/nL, as functions of the temperature. 

The constant Cid corresponds to ideal degenerate nonrelativistic neutron gas at zero temperature. Therefore, at 
T = the second term in the right hand side of Eq. (jXI.2p cancels the term Fid in Eq. (jXI.ip . Thus, at zero temperature 
we have a polytropic EoS (for both the phases) whereas for T > it has more general form defined by Eq. (jXI.f p . 
What why we use the term "quasi-polytropic star" . 

We start hydrodynamic calculations from a hydrostatic configuration at T = 0. According to Eqs. (jXI.f IXI.2IXI.3p 
at zero temperature we have 

FlmH = Cl,h Pl,h = (7l,h - 1) Cl3 n^^-^, (XI.5) 

ML,H = ^^^- (XI.6) 

From Eqs. (|XI.4p it follows that for a polytropic EoS, such as given by Eqs. (|XI.5IXI.6p . the phase transition density 
jump A is unambiguously determined by adiabatic indices of L and H phases: 

^^UH^jHilL-l) ^XI.7) 
TT-L ILKIH - 1) 
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B. Initial model 



Our calculations of the phase transition hydrodynamics were performed in dimensionless variables. Everywhere 
below (in the text and figures) the dimensionless values will be implied if not specified explicitly. 

The following units for the temperature T, density p, velocity V , time t, pressure P, and energy E were used 



T : 



ruu GMs 




■ Ms ' 


, V: 


iGMs 


_ h Rs _ 


, P ■ 




[y Rs \ 



GMs 



P 



'gmY 




\GMn 


47ri?f 


, E: 


[ Rs \ 



(XL. 



where k}, is the Boltzmann constant, G — gravitational constant, Ms and Rs — the mass and initial radius of the 
star, respectively. 

As an example, Mg and Rs for a hot protoneutron star formed in the collapse of the SN core in several tens of ms 
after bounce can be estimated as Ms = 1.4M0 Rs ~ 39 km. Substituting these values for Ms and Rs in Eqs. (jXI.8[) 
we obtain 



[T]«49MeV, [ p ] « 3.7 x 10^^ g/^m^ 



V 



69 000 km/s , 



[ i ] w 0.57 ms , [ P ] w 1.7 X lO'^^ erg/cm^ [ £: ] w 1.3 x 10^^ er, 



(XI.9) 



Instead of the particle number density n, we will use hereafter the mass density p that is connected with n by the 
relation p = rriun where to„ is the atomic mass unit. 

The initial configuration in our calculations was represented by a hydrostatically equilibrium star consisting of 
L- phase matter at T = with the adiabatic index 7^, — 5/3. The dimensionless form of Eq. (jXI.5p for pressure reads 
as 



P = Ap^^^. 



A 



3 



G(47r)2/3My^i?, 



(XI.IO) 



So the initial model is a polytrope of index fc — 1.5 with A — 0.078484 and dimensionless central density pc — 17.97. 
The initial central density at i = was assumed to be just at the beginning of the phase transition. 

In order to simulate the fact that in a core-collapse SN the stellar core is actually rather hot and the pressure can 
exceed that for ideal Fermi gas at zero temperature, the factor Cl was chosen to be equal to 2.83 x Cid- If the density 
jump A and 7^ are specified at T = then one can find 7/f from Eq. (|XI.7p . Next, the factor Ch can be estimated 
from the relation 

PL(0,Pe) =PH(0,Ape). (XI. 11) 

For the main value A ~ 1.525 (at T = 0) in our calculations and 7^ = 5/3 one can derive from Eq. (jXI.7p ~ 1.355. 



In subsequent hydrodynamic calculations the quantities C^^u and 7^^// at any temperatures remain fixed and equal 
to their values estimated at T = as is described above. 

The only parameter that we change in different versions of our calculations is the value of A at T = 0. 

Figure[T7| shows typical dependence of pressure on density at different dimensionless temperatures T = (initial 
model), 1 and 2 (for version A = 1.525 at T = 0). 

An example of the phase diagram in our calculations is shown in Fig.[T8]for A = 1.525 at T = 0. The domains at the 
left of the solid curve and from the right of the dashed one cover the low-density and high-density single-phase states, 
respectively. The mixed-phase states occur everywhere over the region between those curves. Here the pressure does 
not depend on density at given temperature. 

One can see that the density jump A decreases with increasing temperature. 



C. Results of hydrodynamical calculations 



The hydrodynamic calculations were run with fully implicit Lagrangian code, the star being divided onto 1000 mass 
zones. 
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FIG. 17: Pressure as a function of density at different temperatures T = 0, 1, 2 




P 

FIG. 18: The phase diagram for A = 1.525 at T = 

In order to activate the phase transition the hydrostatic equihbrium was perturbed by imposing a small inward 
velocity Vq to stellar matter (Vq = —ar). The total kinetic energy injected in the star by such a perturbation was as 
small as 5 X 10^* of the star gravitational energy Eg for A = 1.525. In this way several models were calculated for 
different values of A at T = 0. Below we discuss the results for the model with supercritical density jump A — 1.525 
at T = 0. 

Figures fT9H2T] show the temporal evolution of density, temperature and velocity from the time of the loss of stability 
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FIG. 19: Temporal evolution of density 



t = (curve labelled 1) on account of the appearance of new phase to the approach to final quasi-equilibrium 
configuration t = 10 (curve labelled 5). Intermediate curves correspond to times t = 6.7(2), t = 7.4(3), and 
t = 7.7 (4). Solid sections of the curves belong to domains contain low-density phase whereas dashed ones correspond 
to high-density phase. All the quantities are shown as the functions of dimensionless mass coordinate < m < 1. 




FIG. 20: Temporal evolution of temperature 



One can clearly observe that the conversion of low-density matter into high-density one has a nature of mini 
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FIG. 21: Temporal evolution of velocity 

collapse and finally results in hydrodynamic transition of the star into new equilibrium. The onset of a new phase 
at stellar center causes pulsational instability. The interface between two phases generates progressive waves that 
excite oscillation of the new-phase core. According to our numerical experiments, the heat-transport processes turn 
out to be of significant importance here for to convert the kinetic energy of pulsations in heat and to redistribute it in 
matter. To include convection in the mixing length approximation (see e.g. |ll7l |) proves to be sufficient in our case. 

As the core of new phase increases the collapse gathers strength and outgoing shock front appears at the phase 
interface. Both the density and temperature are growing quickly behind the front. This process becomes similar to 
standard bounce during the collapse of stellar cores without phase transition. Ultimately the core stops to compress 
and rebounds under the action of ram pressure creating outgoing shock wave that reaches outermost low density 
layers, cumulates and ejects the envelope of a mass Mejc 0.05Ms. 

We have calculated a number of models with different values of the initial density jump A at T = 0. 

Figure[52] shows kinetic energy of ejected matter E^jc m terms of the initial gravitation energy of the star as a 
function of A. For A > |, Eejc has virtually constant value of about 12% of initial gravitational energy since the 
initial configuration is a polytrope of index 1.5 with the dimensionless gravitational energy Eg — —6/7. When A 
decreases below critical value 3/2, E^jc fahs down abruptly — at subcritical A = 1.4 it only amounts to 0.4%. It is 
significant that in the subcritical regime one has to increase the initial inward velocity in order to excite undamped 
(converting finally in collapse) oscillations of the new-phase core. For A = 1.4 the kinetic energy of the initial 
perturbation of about 2 x 10~^|i?g| proved to be sufficient. 

The lower A, the larger the new-phase core should be. This is in full agreement with the work of Blinnikov [ll6| . 
If the mass of a subcritical core is not large enough the core remains stable. Moreover, for A < 3/2 the ejection of 
stellar envelope proceeds owing to numerous pulsations each entailing small leakages of matter rather than a single 
outflow forced by a shock wave. 



D. Discussion 

Our calculations provide an example of hydrodynamical development of the phase transition-induced instability in 
collapsing stellar cores. The onset of the new-phase inner central core gives rise to a shock wave that ejects though low 
mass but of high velocity envelope. For our model the kinetic energy of the ejection can reach about 12% of the initial 
gravitational energy of the core. Taking into account the unit of energy from Eq. ()XI.9p we obtain Eejc ~ 1-3 x lO'^^erg 
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FIG. 22: Kinetic energy of ejected matter versus A at T = 

(A = 1.525). This is even more than necessary to conform to astronomical observations. 

In reahty the final ejection of SN envelope will occur after the following interaction of the phase transition induced 
shock wave with an outer standing accretion shock that decomposes infalling iron-peak nuclei into free nucleons. 
Thus, behind the accreting shock there exists a layer enriched with free neutrons and protons. When crossing this 
layer the phase transition shock has only to slightly accelerate it to cause its expansion accompanied by cooling and 
recombination of neutrons and proton back into heavy nuclei. The energy liberation due to such a recombination 
(about 7 — 8MeV per recombined nucleon) could facilitate the final expulsion of SN envelope [ll8j |. So, the phase 
transition shock can act at least as a trigger of the SN explosion. 

The gravitational energy released in the hydrodynamical transition of the star into a new-phase state converts mostly 
to heat while the density increases by a factor of 10^ (Figs.[T9 l[20l) . Therefore, the processes of heat transport becomes 
of primary importance. The major factors here are expected to be convection and neutrino transport in quark-gluon 
plasma. The latter was neglected in our calculations whereas convection was included. Hence, our explosion energy 
Eejc seems to be underestimated because additional energy supply provided by neutrino heat transport to outer layers 
can only contribute to the ejection energy. 

XII. CONCLUSIONS AND OUTLOOK 

We briefly reviewed in this paper vast material related to the structure of QCD vacuum and the properties of 
neutron stars and presented new results on the physics of the density- induced deconflnement (Sect. Ill), relativistic 
QCB model (Sects. V- VII), EoS of nuclear matter with an admixture of MQS and possible phase transitions (Sects. 
IX-X), and finally and most importantly, we have considered in Section XI the evolution of newly born protoneutron 
stars under the influence of phase transition in nuclear matter. 

We have shown that the phase transition curve in (yit, T) plane of states of the hadron-quark-gluon matter can be 
understood and predicted on the basis of the vacuum energy dominance in good agreement with available lattice 
data, however, the resulting critical density for T < 20 MeV is very high, pc ^ 30po, which leaves room for possible 
intermediate phase transitions in normal nuclear matter below the deconfined quark phase. 

To describe this region, one needs a more adequate description of the nucleon-nucleon forces, which includes quark 
degrees of freedom. We envisaged old QCB model, which was reformulated in modern relativistic form and well 
describes the experimental data. Having at hands new NN (and, in principle, 3iV, 4iV, ..) forces, we have formulated 
new system of equations for the EoS of nuclear matter with an admixture of 6QS, and here the main question is 
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whether MQS play a role of "primitives" , i.e., the instantaneous objects generating hard NN repulsive cores, or they 
exist further in the environment on the mass shell and reveal themselves as a new physical component of high-density 
nuclear matter. We have shown that this new phase may appear when the density increases up to (2 3)po, while 
more quantitative analysis is currently in progress. From this side, one can apply the resulting EoS to describe the 
interiors of neutron stars, the cumulative effect in reactions on nuclei, and study the role of MQS in phase transitions 
in nuclear matter. 

One should emphasize a twofold role of the phase transition of the considered type in physics of supernovae. On 
the one hand, the phase transition could be a missing ingredient necessary to resolve a 30-years puzzle of supernova 
explosion in spherical symmetry. On the other hand, as demonstrated by recent calculations in Ref. [i^, supernova 
outburst induced by the phase transition is accompanied by a second peak on its neutrino light curve that occurs 
in about 250 ms after core-collapse bounce. The temporal and neutrino-flavor properties of such a peak provide 
unique information about the EoS of superdense matter. Thus, nuclear and particle physicists have a good chance to 
experimentally check their findings in case of detection the neutrino signal from stellar core collapse in our Galaxy. As 
for astrophysicists, further detailed and systematic modeling of the phase transition effects in core-collapse supernovae 
is of great importance. 
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XIII. APPENDIX A: THE 6q DESCRIPTION OF ^So AND STATES 



We follow here results of Ref. [49] and the last reference in 5l|. For ^5*0 and '^5'i we neglect Vhh and keep only 
two levels of the QqA4S. The following form of Vhqh is used 

V^^,^(r,r') - [(-gi - E)mir) - c^6{r - - E)r^^{r') - c^5{r' - h)] ^ cl5{r - h)5{r^ - h) ^ ^^^^^ ^^ 

The P-matrix for this potential has the form 

^w-fc+^+#4r+#4r- (XIII.2) 



771 [h) E- El E - E2 
The quantity k is defined due to orthogonality of rj^ (r) to be [5l| 

1 d 



?7i/(r) dr 



k. (XIIL3) 



From the Schrodinger equation for the total wave function ^{r) one obtains /3i 

Then the values of all parameters reproducing the phases with good accuracy in the interval < Tkiu < 515 MeV 
are as follows 

^Sa: b = 7.42 GcV"\ /3i = 0.27 GeV, 

ci = 0.34 GeV^/^ C2 = 0.371 GeV^/^ 
El = 0.23 GeV, E2 = 1 GeV, 
^Si: b ^ 7.025 GeV"\ A = 0.268 GeV, 
ci = 0.343 GeV^/^ ca ^ 0.445 GeV^/^ 
-El 0.243 GeV, E2 = 1 GeV. 

XIV. APPENDIX B: TRANSFORMATION PROPERTIES OF DIBARYON CURRENTS 

In this section, transformation properties of the nucleon wave functions and the operators O entering two-nucleon 
currents under the charge conjugation are given. The nucleon wave functions "i/ = ip <S) x products of bispinors ip 
and isospinors x- 

A. Isospin C-parity 

The isospin C-parity operator and its transformation properties under the main algebraic operations are given by 

Cr = iT^ 



1 

-1 



Cj = CjT^ = C+ = ~Ci 
C*i = Ci. 



The C-conjugation of isospinors is defined by 



Xc - Czx+^, 

xt ^ (Xc)+ = -x'Ci. 
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The isospin hermitian matrices are expanded over the set s = (1,t"), the components of the set transform as 
foUows 



C'l 



Ci = 



1 



(XIV.l) 



The equivalent form reads Cf'E.^Ci = (— where / is isospin of the element. If xi and X2 are g-numbers, 

B. Lorentz C-parity: 



In the standard representation 8J], the C-conjugation matrix and its transformation properties under the main 
algebraic operations are given by 



Cl 



cT 



la^ 
ia^ 



-Cr 



d = Cl=Cl. 
The C-conjugated bispinors are defined as follows 

Matrices acting on the bispinors can be expanded over the set F = (1,75,7'^, 

T 



C 



/I \ 



Cr = 



/I 

757'' 



the components of the set obey 



(XIV.2) 



The bilinear forms are transformed as follows 



the sign can be read off from Eq. (jXIV.2[) . 



C. Combined C-parity 



The C-conjugation operator acting on the nucleon wave functions ~ "ip <Si x lias the following properties: 

C - ClCi, 

= c-^ = c+ = c* = c, 
c = -c. 

The C-conjugated nucleon wave functions are given by 



*e = (*c) = -*'C. 

The matrices O entering two-nucleon currents can be expanded over products of the matrices of S and F. They 
transform according to the rule 



c^s^r^c = ±sr. 
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The bilinear forms are transformed under the permutations as foUows 

The signs are fixed by Eqs. (|XIV.ip and (|XIV.2p . 

D. Even and odd bilinear forms 

Using the above rules, one finds even bilinear forms under permutations of the nucleon fields. These forms are 
constructed using matrices 

(1, il5, 757p). 

The odd bilinear forms are constructed using matrices 

r" (1, i75, 757^) , 

The odd structures including the first-order differential operator, that enter two-nucleon currents, are listed in Table 

13 
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